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OBJECTIVE 
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Study  and  contribute  to  the  advancement  of  integrated  optical  circuits  (IOC’s)  and 
fiber  optics  and  their  application  in  information  transfer  and  processing  systems.  Emphasize 
applications  in  Navy  and  other  DoD  programs  in  airborne,  shipboard,  land-based,  and  under- 
sea systems. 

Contribute  to  and  evaluate  the  state  of  the  art  in  material  and  device  techniques  fo 
the  fabrication  of  IOC  devices  such  as  sources,  couplers,  switches,  and  modulators. 
Determine  the  feasibility  of  employing  fiber  optics  and  IOC  technology  in  military  appli- 
cations such  as  high-capacity  (multi-GHz)  data  links,  data  buses,  and  nonradiating  and 
radiation-resistant  communication  systems. 


RESULTS 

1 . Theoretical  analysis  of  the  propagation  constants  in  diffused  optical  waveguides 
has  been  performed  by  use  of  various  approximating  techniques. 

2.  Eleetrooptic  modulation  of  light  in  a planar  waveguide  fabricated  from  a zinc-blende 
type  crystal  has  been  analyzed. 

3.  Directional  coupling  has  been  observed  in  ZnSe  by  use  of  conventional 
photolithography. 

4.  The  generation  of  closely  spaced  lines  with  extremely  smooth  edges  in  SiOT 
deposited  on  semiconductor  substrates  has  been  demonstrated. 

5.  Characterization  of  Ge-doped  double-heterojunction  light-en.itting  diodes  and 
cw  laser  diodes  has  been  performed. 

6.  Threshold  shifts  in  Ga|.x  A I x As  double-heterojunctior.  lasers  with  a corrugated 
interface  have  been  investigated  theoretically. 

7.  Input/output  waveguide  coupling  structures  have  been  fabricated  and  analyzed. 

8.  Coupling  from  planar  to  linear  waveguides  by  means  ol  horn  structures  has  been 
evaluated  theoretically  along  with  transients  in  periodic  couplers. 

RECOMMENDATIONS 

1 . Develop  optical  interconnections  between  waveguides,  fibers,  and  lasers  using 
horns  or  other  structures. 

2.  Develop  optical  integrated  circuit  elements  such  as  two-port  optical  switches  and 
electrooptieal  modulators. 

3.  Develop  low-loss,  single-mode  glass  fibers,  1-3  km  in  length. 
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4.  Provide  continued  support  to  the  development  of  continuous  room-temperature 
laser  diodes. 

5.  Develop  electronic  drivers  for  electrooptic  switches,  modulators,  and  detector 
amplifiers. 

6.  Assemble  an  IOC  breadboard  demonstration  system  incorporating  the  low-loss 
fibers,  laser  diodes,  and  integrated  optical  circuits  developed  in  this  program. 
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INTRODUCTION 


The  lubrication  of  miniature  solid-state  optical  components  connected  by  channel 
waveguides  is  becoming  feasible  with  the  advancement  ol  such  disciplines  as  the  mater, 
sciences  quantum  electronics,  and  guided-wave  optics.  Integrated  optical  components 
sources  detectors,  modulators,  and  various  coupling  elements  on  one  or  more  tiny  su 
s, rates  will  comprise  systems  much  smaller  in  si/e  and  weight  than  optical  systems  employ- 
ing discrete  components.These  components  will  perform  a number  ol  I unctions  in  the  area 
of  optical  communications.  They  include-  rapid  modulation  and  ^ 

elements  using  applied  fields  to  generate  small  rclractive  index  changes,  coupln  t terntg 
of  s.mials  light  detection  by  p-n  junctions  or  other  structures  in  thin  lilrns.  and  hgh 
geMK'raUon  by  thin-f il m las/r  elements.  The  new  systems  will  be  much  less  susceptible  to 
environ ime ntal  hazards,  such  as  mechanical  vibrations,  extremes  in 
tromagnetie  fluctuations  because  of  their  small  size  and  packaging  density.  In  addition, 
wideband  active  components,  such  as  waveguide  electrooptic  modulators  and  switches, 
:i“o,a£  a,  .cry  low  power  levels  because  ol  ,he  small  dimensions .involved. 

Guided-wave  optical  components  have  found  an  important  potential  use  in  the 
area  of  o pt£? com, min ica, ions  becauw  o.  recen,  progress  in  '''f 
fiber-optic  waveguides  with  very  low  losses,  fiberoptic  wavcguijk-s  wi,  h o w » b * 

I-:  dB/km  a,  0.85-pm  and  I .Ob-nm  wavelengths < or  t.aAs  and  Nd-YA  . C aspec ' yl 
and  single-mode  fibers  having  anticipated  bandw.dths  as  high  as  1 0 GHz  lor  a 1 km , kngin 
immensely  widen  the  horizon  of  optical  communications.  Single-mode  devices  mils  c 
»d  tltnple  energy  efficiently  with  the  libers  and  to  process  the  optical  mlorma- 
.Jon  efficiently  at  rates  approaching  the  bandwidth  capacity  ol  the  libers. 

fiber-optic  waveguide  systems  oiler  significant  advantages  lor  military  inlorma- 

tion  trai'str.  71  nmncdiat,  ly!  with  discrete  components  and 
the  future  with  single-mode  libers  and  integrated  optical  elements.  1 hssc  adsanlag 
inelude  freedom  from  elec, ton, ague, ic  mlerfereucc  If. Ml.  I MP,,  elnmnalion  o grounding 
problems,  and  elimination  of  signal  leakage,  as  well  as  the  potential  lor Targe  » 

sire  weight,  power  consumption,  and  cost.  In  addition  to  high-capauly  point  I o 
eoimiuniications,  a major , mores,  in  inlegralcd  optics  from  a m i, ary  dindpoun  « 
potential  for  implementing  a fiber  optic-lransmisswn  line  mtiltitermmn  (data  b m, 
ilcxing  system  through  low-loss  coupling  and  modulat  ion  elements  This  w ill  l»» 
isolated-terminal,  redundant  informal*,,  transfer.  Urns  lacililahng  the  truly  modular 
(including  distributed  computer)  command  control  and  communications  system, 
rolher  proniising  area  is ,,,  fas,,  high-capacily.  Uigl-dcnsily  nmllipor,  switches  tor 

The  objectives  of  the  program  renorted  here  are  t a advance  the  “^^^N-ivy 

physics  of  integrated  optics  for  military  applications,  to  establish  in  ^ 

and  DoD  programs  a continuing  assessment  ot  system  requirements  . ..  , 

R&D  investments  in  each  application  area,  and  to  produce  prototype  optical  elements 
s^bwlnis  tha"  are  aiincd  " satisfying  these  rerpuremenls.  The  work  on  the  program  lha, 
was  performed  at  Naval  f Ice-Ironies  Laboratory  Center  (NI.LO  and  under  '•““  “'■  ■ 
administered  by  NI.LC  was  in  the  areas  ol  inlegraled-oplical-ciruut  (IOC  I * 1 1 , 

assessment,  materials  for  IOC  devices  and  substrates,  pattern  (abr, canon.  Iheorctiuil  analysis. 

components,  and  system  concepts. 


1 


Ilxpcriincn tally,  effort  continued  to  be  centered  on  the  fabrication  ot  10  components, 
such  as  optical  switches  anil  waveguide  modulators. 

Thi>  report  is  divided  into  the  following  six  sections: 

Materials  and  Fabrication  Techniques  for  Integrated  Optical  Circuitry,  Waveguide 
Device  Properties,  Pattern  Fabrication  Techniques,  Laser  Development,  Fabrication  and 
Analysis  of  Input/Output  Waveguide  Coupling  Structures,  and  Transients  in  Periodic 
Structures:  Coupled  Waves  Approach. 

Within  each  section  are  several  subsections  concerning  specific  work  by  NFLC  or 
others  under  contracts  administered  by  NFLC.  The  performing  organization  is  indicated 
in  the  introduction  to  each  subsection. 
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MATERIALS  AND  FABRICATION  TECHNIQUES  FOR 
INTEGRATED  OPTICAL  CIRCUITRY 


During  recent  years  if  has  become  apparent  that  an  optical  material  analog  to  the 
highly  successful  silicon  substrates  used  in  monolithic  integrated  circuitry  is  presently  not  to 
be  found  in  the  wide  range  of  materials  explored  for  integrated  optical  circuit  (IOC)  use.  Of 
the  materials  investigated1,  the  II -V  compounds,  most  notably  GaAs  and  GaAIAs,  stand  out 
as  the  best  potential  candidates  tor  a totally  monolithic  structure;  that  is,  source,  wave- 
guide components,  and  detector  on  a single  chip,  in  the  medium-to-long-term  development 
time  scale  (5-10  years).  This  technology,  however,  brings  with  it  several  material  and  fab- 
rication problems  which  must  be  solved  before  a complex  IOC  can  be  realized.  Most  notable 
ol  these  is  the  relatively  high  loss  incurred  in  optical  waveguides  at  the  wavelength  of  the 
most  efficient  source.  For  example,  the  leading  contender  as  a source  in  integrated  optics 
is  the  cw  room-temperature  injection  laser  which  emits  at  ~0.85  jum.  The  best  GaAIAs 
waveguides  to  date  have  losses  on  the  order  of  10  30  dB/em  at  that  wavelength.  Choosing  a 
longcr-wavelength  source  such  as  InGaAs  (~1 .0  /am)  with  GaAIAs  as  the  guide  material 
(losses  <5dB/cm  at  1.15  /am)  increases  the  complexity  of  fabrication  and  replaces  the  loss 
problem  with  that  of  lattice  mismatch.  This  is  evident  by  noting  that  the  best  III-V  lasers 
and  waveguides  are  produced  by  planar  liquid-phase  epitaxy  (LPE)  techniques.  Selective 
etching  and  regrowth,  and  multisource  LPF  growth,  however,  do  offer  promise  as  a solution 
to  the  latter  problem.  It  is  apparent,  in  any  event,  that  the  full  benefits  of  monolithic  IOCs 
will  only  be  realized  by  an  extensive  development  of  the  III-V  compounds.  Thus,  it  is  more 
likely  that  integrated  optics  will  follow  an  evolution  similar  to  that  of  the  electronic  IC 
semieondi  '.'tor  technology.  That  is,  the  trend  will  be  towards  a hybrid  approach  in  the 
near  term.  In  the  hybrid  approach,  each  component  group  is  optimized  for  a given 
function,  such  as  light  generation,  active  and  passive  waveguide  components,  or  detection, 
in  terms  of  the  material/fabrication  technique  best  suited  for  the  particular  function.  The 
one  significant  difference  between  electronic  lCs  and  IOCs  is  the  wider  range  of  material 
parameters  required  for  the  latter. 

Table  1 lists  some  of  the  most  promising  material  systems  and  their  associated 
fabrication  methods.  This  is  not  intended  to  be  exhaustive  nor  are  all  possible  fabrication 
techniques  for  the  various  materials  included. 

For  the  near-term  hybrid  technology  of  IOC  development,  diffusion  into  active 
materials  such  as  ZnSe  and  LiNbO^  appears  to  be  the  best  means  of  constructing  active 
waveguides  and  interconnects.  A cw  GaAIAs  injection  laser  is  the  leading  source  for  this 
case,  with  either  InGaAs  or  Si  avalanche  photodiodes  as  the  detectors  in  the  system. 
Optimizing  each  component  in  the  IOC  as  to  performance  places  the  principal  burden  of 
integration  on  coupling  the  various  elements  together.  Until  efficient  high-level  technology 
coupling  schemes  such  as  periodic  or  distributed-gap  couplers  are  reliably  available,  the 
most  likely  coupling  arrangement  will  be  butt-joined  end-fire  coupling.  The  monolithic 
technology  in  IOCs  when  fully  developed  will  of  course  use  continuous  or  distributed 
coupling  between  components. 

To  summarize,  at  this  time  there  is  no  single  material/fabrication  system  which  is 
satisfactory  in  all  respects  for  monolithic  IOC  fabrication.  Considerable  effort  must  be 
devoted  to  material  and  fabrication  techniques  for  a completely  monolithic  technology  to 
be  developed.  Hybrid  systems  optimize  the  individual  components,  but  presently  lack 
optical  coupling  schemes.  For  practical  IOC  devices  a semihybrid  technology  may  prove 
most  efficient  for  ongoing  system  development.  For  the  short  term,  a system  consisting 

1.  See  REFERENCES,  p.  77 
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Material  Representative  Fabrication 


Oi  C.aAlAs*  cw  laser  source,  diffused  waveguides  in  ZnSe,  and  silicon  APD  detection  with 
low-loss,  single-mode  transmission  links  offers  the  highest  overall  capability.  An  ongoing 
task  is  the  assessment  ot  capabilities  offered  by  the  developing  technologies  in  the  entire 
materials/fabrication  area  to  optimize  overall  performance  in  an  IOC  system. 


*Ga AlAs  is  a commonly  accepted  notation  for  GaAl  | xAsx. 
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WAVEGUIDE  DEVICE  PROPERTIES 


PROPAGATION  CONSTANTS  IN  DIFFUSED  OPTICAL  WAVEGUIDES 

I h is  section  reviews  some  theoretical  work  performed  at  NELC  on  the  determination 
ol  waveguide  properties  in  a guide  containing  a graded-index  boundary.  The  exact  propa- 
gation constants  obtainable  in  a guide  containing  an  index  gradient  such  as  that  produced 
by  dillusion  are  essential  lor  optimum  device  performance.  Device  design  parameters  are 
dependent,  sometimes  critically  so.  on  the  exact  nature  of  the  propagation  of  light  in 
the  waveguide.  I he  determination  of  propagation  constants  in  a planar  graded-index 
waveguide  is  a necessary  lirst  step  to  fully  understanding  propagation  in  more-complex 
structures. 


LLECTROOPTIC  MODULATION  OP  LIGHT  IN  A PLANAR  WAVEGUIDE  FABRICATED 
FROM  A ZINC-BLENDE  TYPE  CRYSTAL 

Previous  work  on  waveguide  eleetrooptie  modulators-  at  NT.  LC  dealt  with  the 
fabrication  ol'  devices  in  I l-Vl  compounds  such  as  ZnSe.  vv  Inch  has  the  cubic 
/inc-blende  tv  pc  structure.  This  section  presents  theoretical  calculations  performed  at 
NELC  on  certain  aspects  of  device  performance  which  arise  from  the  confinement  of  the 
optical  energy  in  a waveguide  compared  to  bulk  crystal  electrooptic  modulation. 

Planar  optical  waveguides  have  been  fabricated  in  the  ll-VI  compounds  by  dif- 
lusion-V  I he  dill  used  planar  optical  waveguides  are  fabricated  by  diffusing  an  appropriate 
element  into  single-crystal  ll-VI  compound  semiconductors  to  produce  a composition 
gradient  at  the  surface  of  the  crystal.  In  those  mixed  crystals  in  which  the  refractive  index 
is  a function  of  composition,  a surface  layer  of  sufficiently  high  refractive  index  to  guide 
light  is  formed.  1 he  fabrication  of  optical  waveguides  by  this  technique  has  been  demon- 
strated for  the  diffusion  of  Se  into  CdS,  ( d into  ZnSe,  CM  into  XnS,  and  Se  into  ZnS^. 

II  these  waveguides  are  to  be  employed  in  IOC's,  a knowledge  of  the  propagation 
velocities  of  the  modes  supported  by  these  structures  is  necessary.  The  refractive  index 
prolile  lor  these  waveguides  has  been  measured-'  and  found  to  be  the  error  function  com- 
plement in  agreement  with  simple  diffusion  theory.  The  fields  in  such  waveguides  and  the 
transcendental  equation  for  the  propagation  constants  cannot  be  formulated  in  closed-form 
expressions  involving  only  simple  functions.  Approximate  values  for  the  propagation 
constant  can  be  obtained  by  WKB  techniques  or  from  the  exact  solutions  for  other  refrac- 
tive index  profiles.  Closed-form  solutions  exist  for  square*',  linear7,  anil  exponential^ 
refractive  index  proxies. 

I wo  diffusion  profiles  are  of  interest  for  diffused  waveguides.  The  diffusion  ol  a 
dopant  from  a large  volume  of  gas  into  a crystal  is  described  by  the  model  of  diffusion 
from  a source  of  constant  concentration  into  a semi-infinite  body.  The  profile  of  the 
diflusant  concentration  for  this  ease  is  an  error  function  complement.  If  the  diffusion  is 
Irom  an  minutely  thin  layer  of  the  diflusant  that  is.  an  exhaustible  source  the  con- 
centration profile  of  the  diflusant  is  Gaussian1^.  The  refractive  index  profiles  of  interest 
are  illustrated  in  figure  I.  Each  piofile  can  be  expressed  as 

N(X)  = Ns  + AN  X l'(X). 
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Figure  I . Refractive  index  vs  depth. 

where  Ns  is  the  refractive  index  of  the  substrate  and  Ns  + AN  the  refractive  index  at  the 
surface.  The  function  f(x)  is  given  in  figure  1 in  association  with  the  corresponding  retractive 

index  profile.  . . , , , , . , , . 

Numerical  solutions  for  the  modes  propagating  in  the  waveguides  were  obtained  by 

dividing  the  planar  structure  into  three  regions:  I , the  region  outside  the  substrate;  2,  the 
region  in  the  substrate  close  to  the  surface  where  there  is  significant  spatial  variation  in  the 
refractive  index;  and  3,  the  bulk  of  the  substrate.  For  Tb  modes  the  electric  fields  in  regions 
I and  3 are  given  by  the  expressions 

E|  = Aj  exp  (gj X + (w  t -0  Z| 
and 

i‘3  = A3  exp  |-g3X  + (w  t - 0 Z| , 


where 


o'*  1 T 

gr  = 0“-k“ 


kr  = 4*  nr 


r=  1.3 


and  Ai , Aj  arc  constants.  The  electric  field  in  region  2 is  given  by  the  solution  of  the  differ- 
• in 

ential  equation  lu 

1 1 ( X ) = l k-(X)  -0-|  H ( X ) = 0 
k(X)=  -y-  N(X) 

Hie  magnetic  field  along  the  direction  of  propagation  in  ali  three  regions  is  given  by  the 
expression 


H(X)  = 


h.'(X) 


Solutions  were  obtained  by  assuming  a value  tor  B and  lor  the  magnitude  of  the 
electric  field  in  region  3.  The  continuity  oi  the  electric  and  magnetic  tields  across  the  2,  3 
boundary  was  used  to  obtain  the  magnitude  of  the  electric  field  and  its  tirst  derivative  in 
region  2 at  the  2,  3 boundary.  The  ordinary  differential  equation  for  the  electric  field  was 
integrated  across  region  2 to  the  1 , 2 boundary  by  use  of  a Miln  predictor-corrector  algorithm. 
The  electric  and  magnetic  fields  at  this  boundary  in  regions  1 and  2 were  then  compared.  A 
two-point  interpolation  iteration  technique  was  used  to  find  the  values  ol  g that  matched 
the  fields  at  the  1 . 2 boundary.  The  boundary  2.  3 was  placed  at  the  value  of  X for  which 
f(X)  equals  0.1  X f(o). 

The  accuracy  of  the  computation  was  cheeked  by  comparing  the  propagation  con- 
stant computed  for  a rectangular  profile,  by  use  of  the  above  algorithm,  with  the  exact 
solution.  The  values  of  g agreed  to  three  signit icant  places  and  ot  B to  five  places. 

The  lateral  confinement  of  a mode  in  a dielectric  waveguide  is  analogous  to  the 
quantum  mechanical  model  of  a particle  in  a potential  well.  The  WKB  approximation  has 
been  extensively  developed  for  this  later  problem1 1 . The  technique  is  directly  applicable  to 
the  present  problem  and  yields 

I \/k -<X)-0“  dx  = (n-l/2 in 


as  the  determining  equation  for  the  propagation  constant  0 where  n - 1.2...  and  7 is  the 

point  at  which  the  magnitude  of  |k(X)-  -(J- 1 equals  zero. 

The  propagation  constants  for  the  lowest-order  mode  of  waveguides  with  Gaussian, 
exponential,  square,  and  linear  refractive  index  profiles  are  plotted  in  figures  2 and  3.  where 

«=  ^-(f) 

(-f h 

In  figure  2.  AN  = \'/<  Ns  for  all  waveguides.  In  figure  3.  the  AN  for  each  profile  has  been 
adjusted  in  a manner  such  that  the  integrated  excess  retractive  index. 

/ |N(X)-NS]  dX. 

equals  0.01  X N1  X D for  all  waveguides  where  N1  - AN  for  the  square  profile. 


8 


NORMALIZED  PROPAGATION  CONSTANT  (B)  NORMALIZED  PROPAGATION  CONSTANT  (B) 


Figure  2.  Propagation  constant  for  square,  Gaussian,  linear,  exponential,  and 
error  function  complement  refractive  index  profiles,  AN  = 19?  Ns- 


Figure  3.  Propagation  constant  for  square,  Gaussian,  linear,  exponential,  and 
error  function  complement  refractive  index  profiles,  integrated  refractive 
index  = 0.01  X Ns  X D. 
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In  the  case  ot  the  symmetric  square  waveguide,  the  differential  equations  for  the 
electric  and  magnetic  fields  can  be  normalized  by  using  the  reduced  quantities 

h=  - Ns  K- 

(Ns  + AN)-  - K- 

and 

(I  = |<Ns  + AN)--N“|I/:  I)/X 

Iliese  quantities  do  not  normalize  the  corresponding  equations  when  the  refractive  index 
profile  is  not  constant  with  depth  in  the  guiding  region.  Nevertheless,  calculations  showed 
that  for  the  range  of  Ns  from  2.3  to  3.0  and  AN  from  0.17  to  107  Ns  the  use  of  the 
above  reduced  quantities  did  bring  the  curves  for  B versus  D into  coincidence  to  within  two 
significant  places.  The  propagation  constants  for  waveguides  fabricated  in  ZnTe,  CdS,  CdSe 
and  ZnS  and  measured  at  wavelengths  other  than  6328A  can  therefore  be  found  from 
figures  2 and  3 by  using  the  expressions  for  b and  d. 

It  can  be  seen  from  figures  2 and  3 that  each  of  the  profiles  could  be  used  as  an 
approximation  for  estimating  the  magnitude  of  the  propagation  constant  of  another  of 
the  proliles:  however,  such  an  estimate  would  be  only  a rough  approximation,  The  error 
function  complement  and  the  exponential  profiles,  compared  on  the  basis  of  integrated 
retractive  index  profiles,  yielded  closely  similar  values  of  the  propagation  constant  far 
Irom  cut-off.  At  cut-off  tile  discrepancy  is  large,  and  the  exponential  approximation  yields 
poor  estimates  foi  the  wave  dimensions  at  cut-off. 

The  propagation  constant  for  the  error  function  complement  profile  was  estimated 
by  the  WKB  approximation  and  is  included  in  figures  2 and  3.  It  is  simpler  to  estimate  the 
propagation  constant  by  this  technique  than  by  any  of  the  exact  solutions,  and  the  accuracy 
should  be  adequate  for  many  applications. 

Two  realizations  o'-  a planar  waveguide  modulator  are  illustrated  in  figure  4.  The 
voltage  applied  between  the  electrodes  induces  an  electric  field  predominantly  along  the 
x-axis  in  the  transverse  electrode  configuration^)  and  predominantly  along  the  y-axis 
in  the  parallel  electrode  configuration  (4b).  While  it  is  possible  to  apply  the  electric  field 
along  the  z-axis,  the  significant  increase  in  the  applied  voltage  required  for  this  configura- 
tion eliminates  it  from  practical  consideration. 

The  Planar  waveguide  modulators  discussed  here  are  intrinsically  phase  modulators. 

If  no  static  electric  field  is  applied  to  the  crystal,  then  TK  and  TM  modes  propagate  along 
the  planar  waveguide.  If  a voltage  is  applied  to  the  electrodes,  the  resulting  electric  field 
perturbs  the  dielectric  tensor.  The  modes  then  may  become  mixed  and  the  propagation 
constants  of  some  or  all  of  the  modes  present  may  change  as  a function  of  the  applied 
field.  Phase  modulation  will  result  if  a mode  whose  propagation  velocity  is  sensitive  to  the 
modulating  voltage  is  launched  into  the  waveguide.  Amplitude  modulation  is  achieved  by 
launching  two  modes  with  equal  initial  phase,  usually  a TP  and  a TM  mode,  into  the  modu- 
lator. If  the  sensitivities  of  the  propagation  velocities  of  these  two  modes  to  the  modulation 
voltage  are  different,  then  the  interference  of  the  two  modes  at  the  exit  of  the  modulator 
will  yield  a signal  which  is  amplitude  modulated  by  the  electric  field  applied  to  the  electrodes. 

These  modulators  are  the  waveguide  analog  of  the  well  known  bulk  modulators. 

A compensator  must  be  added  to  the  output  of  the  modulator  to  achieve  maximum 
amplitude  modulation.  The  compensator  can  be  eliminated  if  the  length  of  the  modulator 
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IS  selected  to  be  such  that  the  two  modes  are  in  phase  at  the  exit  of  the  modulator  when  the 
voltage  aeross  the  electrodes  is  zero. 

An  idealization  ol  the  geometry  of  the  modulators  is  shown  in  figure  4c,  where  the 
propagation  is  assumed  to  be  in  the  z direction.  An  eleetrooptically  active  layer,  2d  thick, 
is  bounded  above  and  below  by  nonactive  regions  of  lower  refractive  index.  It  is  assumed 
in  addition  that  there  is  no  variation  in  the  properties  of  the  modulator  or  in  the  electro- 
magnetic fields  along  the  y-axis.  This  latter  assumption  would  not  be  a good  idealization 
ol  a single  transverse  mode  channel  waveguide  modulator.  However,  the  idealization  is 
useful  for  those  modulators  described  in  reference  12  which  have  much  greater  width  than 
depth  and  which  support  a very  large  number  of  transverse  modes. 

Under  suitable  circumstances  the  crystallographic  axis  could  have  an  arbitrary 
oiicntation  with  respect  to  the  waveguide  axis.  Any  given  fabrication  technique,  however, 
usually  yields  only  a limited  number  of  possibilities.  For  the  cubic  crystals  the  principal 
axes  in  the  absence  of  an  applied  field  are  aligned  with  the  crystallographic  axes,  and  the 
index  ellipsoid  is 


where  double  prime  represents  the  crystallographic  axes. 

In  the  presence  of  an  applied  electric  field  £x,  £y,  £z,  the  index  ellipsoid  deforms  to 
become*’ 


x»2  + y"2  + /»2 


-> 

11“ 


~ + 2'4|  (&xy''z"  + £yz"x”  + £zN"y")  = 1. 


Hie  axes  ol  the  deformed  ellipsoid,  the  primed  axes  of  figure  5,  arc  rotated  from  the 
crystallographic  axes.  For  purposes  of  illustration  the  first  two  cases  examined  by  Namba^ 
are  illustrated  in  figure  5. 

The  dielectric  tensor  in  the  primed  coordinate  system  is 


Rotation  of  this  tensor  to  the  unprimed  or  waveguide  axis  eoordinate  system  yields  the 
tensor 
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H ep  and  03  equal  zero,  the  modes  that  can  propagate  will  remain  TF  and  TM  modes. 
The  simplest  example  of  this  case  occurs  if  the  waveguide  axis  coincides  with  the  principal 
axis  of  the  deformed  ellipsoid.  Since  the  zinc-blende  crystals  cleave  along  the  (1 10)  planes, 
a waveguide  oriented  as  in  figure  6 is  easily  labricated.  Such  a waveguide  used  in  a trans- 
verse modulator  becomes  ease  1 ot  tigure  5.  The  dielectric  tensor  in  the  unprimed  coordi- 
nate system  for  this  orientation  of  the  crystal  is  1 3. 

(n  + p)- 

(n  -p)~ 

where 

p = '/in  ^ r^  j 6 


From  Maxwell’s  equations  the  differential  equations  for  Fy  and  Hy  within  the  guiding 
region  are 


Kv"  + (k^  - B-)Fy  - 0 , 


and 


Hy"+ej  (k|  - B-)  Hy  = 0, 


(1) 

(2) 


Figure  6.  Crystal  and  optical  waveguide  orientation. 
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where 

k i 2 = k22  = and  B is  the  propagation  velocity. 

These  equations  are  uncoupled  and  the  modes  remain  TE  and  TM. 


For  the  lowest-order  TE  mode  let  Ey  be 
Ey  = A cos  Ux 

within  the  guiding  region  and 

Ey  = A cos  Ude*v(x-d) 

in  the  evanescent  region.  The  characteristic  equation  for  this  mode  is 
U sin  Ud  - V cos  Ud  = 0, 


(3) 


where 

U2  + v2  = K22  = k22-k02, 

k02  = uVq, 

and  cq  >s  the  dielectric  constant  of  the  bounding  regions. 

The  corresponding  characteristic  equation  for  the  lowest-order  TM  mode  is 

e2U  sin  Ud  - eg  8 cos  Ud  = 0,  (4) 


where 


f3  U2  + v2  = K2  = k2-k2 
C1 

The  above  equations  were  solved  numerically  for  the  propagation  velocities  ot  the 
lowest-order  TE  and  TM  modes  as  a function  of  the  static  electric  lield.  The  waveguide 
used  in  the  numerical  calculations  had  the  properties 

d = 0.25  pm, 
n0  = 2.58, 
n = 1.01  X N0, 

T4i  = 2 X 10"'0  cm/volt, 

and 

X = 632.8  nm. 
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Tli  esc  values  were  chosen  in  order  to  model  the  optical  waveguides  fabricated  in  ZnSe^  * 
and  tested  with  a HeNe  laser,  The  magnitude  of  d is  not  critical  for  this  purpose. 

The  variation  of  the  propagation  velocity  of  the  TE  mode  with  the  magnitude  of 
the  static  electric  field  is  plotted  in  figure  7,  It  is  evident  from  this  figure  that  the  variation 
of  the  propagation  velocity  of  this  mode  is  closely  linear. 


& (STATIC  ELECTRIC  FIELD.  105  VOLTS/cm) 
figure  7.  Differential  phase  delay  for  Tl£  and  TM  modes,  case  I. 


The  derivative  of  the  propagation  velocity  with  respect  to  the  applied  electric  field  is  given 
by  the  expression 

, 2 

dB  V-  + VK|d  dk|  dkj 

dfi"F  Ky  + VK^d  di"  _ ' (IF 


where  the  factor  is  the  rate  of  change  of  the  propagation  velocity  for  a plane  wave  in 
diu 


the  bulk  material  propagating  in  the  same  direction  as  the  waveguide.  The  factor  f varies 
from  zero  to  unity  as  v varies  from  its  lower  bound  of  zero  to  its  upper  bound  of  K j . The 
change  in  the  velocity  of  propagation  in  the  waveguide  is  then  always  less  than  the  change  in 
the  propagation  velocity  in  the  bulk  material. 


The  variation  of  the  propagation  velocity  of  the  TM  mode  with  applied  field  is 
also  shown  in  figure  7.  The  propagation  velocity  of  the  TM  plane  wave  in  the  bulk  material 
oriented  in  the  same  direction  as  the  TM  mode  in  the  waveguide  is  independent  of  the 
applit  I voltage. 

As  the  second  example,  consider  the  electric  field  to  be  applied  transversely  for  the 
same  waveguide  orientation  as  before.  The  dielectric  tensor  is 


in  the  primed  coordinate  system  and 


in  the  unprimed  coordinate  system.  The  equations  for  Ey  and  Hy  for  this  case  are 

F.y"  + j^k2  - B2  - ^(2pn)2  J Ey  = gjB*i  Hy  (5) 

and 

Hv"  +—  (k2  - B2)  H„  = -ojB^  (2pn)  Ey. 
y C]  j e | y 

where 

k2  - co 2 pi  (n2  + p2). 

These  equations  are  coupled,  and  the  modes  that  propagate  are  no  longer  pure  TE  and  TM 
modes  but  rather  mixed  modes. 

These  mixed  modes  are  not  coupled  but  are  each  orthogonal  to  all  other  modes  that 
may  propagate.  The  departure  from  being  pure  TE  and  TM  modes  is  very  slight  if  the  pro- 
pagation velocities  of  the  TE  and  TM  modes  of  the  same  order  are  not  closely  equal.  In  the 
absence  of  a voltage  applied  to  the  electrodes,  there  is  no  coupling,  and  the  modes  are  TE 
and  TM.  Only  with  the  application  of  a voltage  do  the  modes  change  to  what  will  be  termed 
here  QTE  and  QTM  modes  for  quasi  TE  and  quasi  TM.  The  QTE  mode  then  grows  from  the 
TE  mode  and  the  QTM  mode  from  the  TM  mode  with  the  application  of  a voltage  to  the 
electrodes. 

The  changes  in  the  propagation  velocities  of  the  QTE  and  QTM  modes  are  plotted 
versus  the  transverse  electric  field  intensity  in  figure  8.  The  zero  field  intensity  propagation 
velocity  difference  has  been  suppressed.  These  numerical  results  are  for  an  optical  wave- 
guide of  the  same  properties  as  the  previous  numerical  example.  The  variation  of  the 
propagation  velocity  with  applied  field  is  square  law  to  within  5'/( . The  change  in  the  pro- 
pagation velocity  with  applied  voltage  in  the  numerical  example  can  be  approximated  by 
the  expression 

BQTE(i:,-BQTM(S,  = (BTE-BTM>=  l-48X  lO"13^2.  <6> 
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CASE  2 


(STATIC  ELECTRIC  FIELD.  104VOLTS/cm> 

™5' OTH  “US'S  t"  **  B C,C'"ic  “d  fl"  0TE 


Th 


up p roximate  expression 


BQT1  ( ) -BqTM(G  ) = (BT].  -BTM) 


4 

k\  BT1.n  r4f 
wUc(Uni  - Ue) 


fi2 


(7) 


lor  the  propagation  velocity  difference  is  derived  in  reference  15.  For  the  waveguide 

w in  mT  3|’PrOXi",a,'°"  r°r,hc  covl  lic'tfnt  or  the  6 2 «.-<Ilia|s  3.74  x 10-13 

I n n,  r°rir'1'  “ r ,i<,n  J<K'S  ,t|,,odu“  *"«  -l^dr.Uc  dependence  upon  8.  ' 
k magnitude  ol  the  coefficient  is  in  error  by  a factor  of  2.5. 

In  this  second  case  the  results  for  the  waveguide  modulator  are  substantially 

d‘  erent  than  lor  the  hulk  modulator.  The  variation  of  the  propagation  velocity  difference 

j app  ic  vo  tage  is  quadratic  in  the  waveguide  modulator  and  linear  for  the  bulk 

I he  orientation  of  the  polarizers  used  for  amplitude  modulation  is  also  different 
i ie  waveguide  modulator  the  waveguide  boundaries  dominate  over  the  weak  induced 
anisotropy  and  the  boundaries  define  the  T1  and  TM  modes.  The  polarizer  oriental^ 
i amphu.de  modulation  is  45°  to  the  waveguide  normal  or  the  wave- 

XU,cd  by  u s,snal  P°lar,zcd  ^ <1 1 1>  or  <TT1>  direction.  In  the  bulk 


18 


modulator  the  weak  anisotropy  defines  the  or  entation  of  the  TE  and  the  TM  modes 
in  a manner  such  that  the  polarization  of  the  incident  light  must  be  along  the  <0(H>or<l  I0> 
direction  to  achieve  maximum  amplitude  modulation.  This  direction  is  then  at  45  to 
that  required  for  the  waveguide  modulator.  In  contrast  to  these  results  the  fields  of  the 
QTL  and  QTM  modes  deviate  only  slightly  from  being  TE  and  TM  modes,  respectively. 

The  magnitude  of  the  light  intensity  out  of  the  modulator  is  related  to  the  differ- 
ential phase  delay  by  the  expression 


I = cos 


2 


ns) 

2 , 


where  the  differential  phase  lelay  is 


r(£)  = [Bqte(£)-Bqtm(£)]  L (8) 

at  wavelength  X for  a modulator  of  length  L. 

Experimental  results  verity  the  polarizer  orientation  and,  at  low  modulation  voltages, 
the  quadratic  dependence  (11).  Experimental  data  for  the  amplitude  of  the  modulator 
output  versus  the  applied  voltage  are  reproduced  in  figure  9.  For  comparison  the  magnitudes 
of  the  light  outputs  for  a linear  and  a quadratic  propagation  velocity  difference  dependence 
on  applied  field  are  also  shown. 
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figure  9.  Modulator  output  vs  electrode  voltage  (relative  units),  measured  data, 
linear  phase  dependence,  and  square  law  phase  dependence. 
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Numerous  other  special  cases  ol  the  planar  waveguide  modulator  could  he  considered. 
However,  the  two  cases  discussed  here  illustrate  the  new  phenomena  that  can  be  introduced 
by  the  presence  ot  the  guiding  structure.  In  addition,  the  waveguide  modulators  that  are 
projected  to  be  useful  in  integrated  optical  circuits  are  single-mode  channel  waveguides. 

Hie  strong  lateral  confinement  ol  such  waveguides  greatly  complicates  the  analysis,  and 
the  approximate  analysis  which  yields  such  simple  expressions  for  the  case  considered  here 
is  not  applicable. 


EXPERIMENTAL  FABRICATION  OF  FREQUENCY-SELECTIVE  WAVEGUIDES  IN 
GLASS 

Previous  work  ^ " at  Nl  l.(  on  the  lubrication  ol  optical  waveguides  in  glass  by 
use  ol  a focused  cvv  laser  beam  has  resulted  in  low  loss  waveguides.  Flic  process  involves 
c\v  laser  heating  of  a glass  that  is  highly  absorbing  at  (he  laser  wavelength  but  has  very 
low  losses  at  the  longer  wavelengths  used  tor  waveguiding.  This  section  reports  progress 
in  waveguide  fabrication  by  cvv  laser  heating. 

Attempts  were  made  to  demonstrate  Irequency-scicctivc  waveguiding  in  channel 
optical  waveguides.  A periodic  slructure  was  imposed  on  the  waveguide  by  chopping  the 
ew  laser  beam  by  means  of  a mechanical  chopper. 

Distributed  feedback  should  be  observable  in  waveguides  possessing  a period  p 

p=3iT  (2N  + n 

with  N - 0,  1,2,...  The  fabricated  waveguides  were  about  3 pm  wide.  Orders  with  N = 0, 
1 , 2 were  fabricated,  employing  different  laser  powers,  chopping  frequencies,  and  ex- 
posure times. 

To  increase  the  sensitivity  of  the  detection  method,  a Princeton  Applied  Research 
lock-in  amplifier  in  the  differential  mode  was  used.  The  chopped  light  output  from  the 
near  field  guiding  area  was  collected  on  a silicon  diode,  giving  rise  to  signal  A.  Scattered 
light  in  the  vicinity  ot  the  waveguide  was  collected  on  another  silicon  diode,  resulting  in 
signal  B.  With  the  aid  of  a variable  transmission  filter  the  amplitudes  of  signals  A and  B 
were  adjusted  until 

At X ) - B(X)  ~0. 


The  difference  signal  was  then  amplified.  At  p - Xj,  one  would  expect  strong  reflection  of 
this  wavelength  due  to  the  corrugated  structure  of  the  waveguide,  resulting  in  a signal 

A(X|)-B(X|)#0. 

In  performing  these  experiments,  it  became  apparent  that  a very  high  level  of  noise 
was  obstructing  any  observable  effect. 

C onsiderable  effort  was  spent  to  reduce  noise,  which  had  its  main  sources  in  many 
mechanical  vibrations  within  the  system.  New  holders  for  mirrors,  objectives,  mountings, 
etc.  had  to  be  specially  constructed  in  order  to  eliminate  the  mechanical  instabilities. 


To  date,  in  spite  of  significant  improvements  in  the  mechanical  and  electronic 
stability  of  the  measuring  apparatus,  no  frequency-dependent  transmission  has  been  observed. 
It  is  felt  that  the  heating/sublimation  mechanism  by  which  the  waveguide*  are  formed  may 
make  small  dimensional  changes  very  difficult  to  accomplish  due  to  the  thermal  con- 
ductivity of  the  glass. 

Work  on  frequency  selectivity  is  now  receiving  lower  priority  due  to  the  stert-up 
of  work  on  fabricating  inout/output  coupling  structures  by  use  of  cw  laser  heating 
techniques. 


PATTERN  FABRICATION  TECHNIQUES 


Fabricating  good-quality  optical  waveguide  components  by  diffusion  requires 
masking  techniques  which  have  extremely  good  pattern  definition  and  repeatability.  At 
present,  two  such  techniques  are  being  investigated  in  this  program  for  this  purpose.  NFLC 
is  purs  ling  state-of-the-art  photolithography  for  producing  the  diffusion  masks  in  SiO->  and 
Hughes  Research  Laboraioiies  is  investigating  electron  beam  lithography.  The  following 
two  sections  describe  the  current  status  of  the  two  approaches 


PHOTOLITHOGRAPHIC  FABRICATION  OF  COUPLED  WAVEGUIDES 

To  date,  photolithography  appears  to  be  the  most  useful  masking  technique  for  con- 
structing IOC  components  by  diffusion  in  the  1I-VI  compound  semiconductors.  Recent  work 
has  resulted  in  significant  improvement  in  the  quality  level  of  simple  waveguide  components 
and  much  promise  with  more-complex  structures.  Basically,  contact  photolithography 
requires  very  smooth  surfaces  for  optimum  reproduction  of  the  photomask  pattern.  Thus, 
large-area  (>lcm~)  substrates  with  extremely  smooth  ('vX/SO  rms)  and  flat  (,wX/2)  surfaces 
are  required  lor  the  best  possible  pattern  resolution.  ZnSe  used  at  NFLC  for  most  of  the  recent 
dill  used  waveguide  work  has  a strong  tendency  to  twin  along  (111)  planes.  The  orientation 
ol  the  substrates  used  in  the  work  has  been<IOO>.  This  results  in  the  intersection  of  the  (111) 
twin  planes  with  the  (100)  top  surface,  resulting  in  very  irregular  surfaces  after  chemical 
polishing  due  to  etch  rate  variations  in  the  vicinity  of  the  twinned  areas.  The  size  of  the 
untwinned  areas  on  a (100)  crystal  was  typically  3 by  10  mm,  although  10  by  10-mm  areas 
could  be  obtained  on  rare  occasions. 

Discussions  with  Y.  S.  Park  and  coworkers  at  Wright-Patterson  Air  Force  Base 
Solid  State  Physics  Laboratory  resulted  in  the  change  to  (1 1 1)  oriented  crystals.  In  this  orien- 
tation the  face  of  the  substrate  is  parallel  to  the  twin  planes,  resulting  in  defect-free  crystals 
ot  20  by  20- mm  area  with  an  average  size  of  10-15  mm  square.  This  is  a significant  increase 
in  usable  area  (~300^)  and  has  resulted  in  rapid  progress  in  photolithographic  pattern 
fabrication. 

Arrays  of  optical  waveguides  of  2-pm  width  separated  by  2-gm  spaces  over  1 cm  long 
have  been  Lbi'cated  in  ZnSe  by  the  diffusion  of  cadmium.  Figure  10  shows  oscillograph 
traces  of  the  light  intensity  at  6328  A detected  by  scanning  a silicon  diode  across  the  output 
ot  an  array  of  4-gm,  single-mode  coupled  waveguides  separated  by  5-pm  spaces.  The  wave- 
guides are  coupled  by  their  evanescent  fields. 

The  intensity  distribution  among  two  coupled  waveguides  can  be  expressed  as 

1 1 (/.)  = cos  - K/. 

h (/.)  = sin  - Kz 

where  / is  the  propagation  length  and  K is  the  coupling  coefficient.  In  an  infinite  array  of 
coupled  waveguides  the  light  intensity  in  the  nth  guide  away  from  the  initially  excited 
(zero1*1)  guide  is  (I  5) 

l„  (z)  = Jn2  (2Kz.)  . 

1 igure  I I shows  the  relative  intensities  of  the  first  few  orders.  From  figure  1 0 the  value  of  K 
is  calculated  to  be  approximately  0.5  mm'1  at  6328  A.  In  addition,  no  scattering  is  visible 


Figure  10.  Directional  coupling  between  parallel  waveguides.  Waveguide 
dimensions  4-//m  width  by  3/am  spaces.  Sample  length:  (a)  1.70  mm. 

(b),  2.70  mm,  (c).  3.55  mm.  Waveguides  fabricated  by  cadmium  diffusion 
in  ZnSe.  Diffusion  depth  2 gm. 


Figure  1 1 . Magnitude  of  the  first  three  orders  of  the  ordinary  Bessel  function 
as  a function  of  argument  K/. 
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from  these  waveguides,  indicating  scattering  losses  of  less  than  1 dB/cm,  which  should  be 
compared  to  the  more  than  20  dB/cm  of  the  previously  reported  directional  couplers  in 
(la As  (16)  fabricated  by  ion  implantation. 

The  significance  of  this  work  is  not  just  that  directional  couplers  were  fabricated,  but 
that  fabrication  was  accomplished  by  conventional  photolithographic  techniques  which  are 
inexpensive  and  quite  straightforward  to  implement.  More-complex  active  structures  arc 
imminent  and  await  only  the  acquisition  of  suitable  masks. 

ELECTRON  BEAM  LITHOGRAPHY 

Hughes  Research  Laboratories  (URL)  has  developed  several  novel  techniques  for  the 
fabrication  of  patterns  suitable  for  the  diffusion  fabrication  of  optical  waveguides.  This  work 
is  oriented  toward  the  generation  of  oxide  diffusion  masks  on  semiconductor  substrates  and 
employs  a computer-controlled  scanning  electron  microscope.  Recent  work  at  HRLhas  con- 
centrated on  the  generation  of  closely  spaced  lines  with  extremely  smooth  edges  etched  in  a 
SiCH  masking  layer  deposited  on  a semiconductor  substrate. 

Two  groups  of  four  samples  each  of  ZnSe,  oriented  ( 1 00),  each  with  a 1 000-A  layer 
of  sputter-deposited  SiCH.  were  supplied  by  NELC  for  electron  beam  exposure  and  pattern 
delineation  in  the  oxide  layers  that  are  used  subsequently  as  diffusion  masks  for  creating 
optical  waveguides.  The  surfaces  of  some  of  these  samples  were  complicated  by  the  presence 
of  several  large  crystals  terminating  in  the  surface.  Patterns  were  delineated  as  much  as  possible 
in  smooth  areas  between  crystal  surface  structures,  which  were  usually  confined  to  a particular 
crystallite.  However,  in  some  cases  the  patterns  intercepted  such  structure. 

Figure  12  shows  a cross  section  of  the  patterns  or  lines  fabricated  with  the  character- 
istic dinvMisions  A (line  width)  and  B (line  spacing).  Chemical  etching  of  the  SiOp  was  used 
after  electron  beam  exposure  and  development  of  the  electron  resist.  This  decision  was  based 
on  successful  results  in  the  past  and  the  fact  that  the  particular  line  widths  required  for  this 
program  are  relatively  large  (>1  pm).  The  validity  of  the  choice  is  shown  in  figure  13,  which 
is  a scanning  electron  micrograph  (2000X)  of  a typical  set  of  three  etch  channels  in  the  SiCb 
films.  Note  tiiat  the  edge  roughness  is  undetectable,  even  at  this  magnification  (probably 
<500  A)- 


Figure  12.  Cross  section  of  the  line  sets  fabricated,  illustrating  the 
two  characteristic  dimensions  A and  13. 
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Figure  13.  Scanning  electron  micrograph  (2000X) 
of  chemically  etched  channels  in  SiOi.  (Channel 
width  is  4.5  /am.  Edge  smoothness  is  better  than 
500  A except  for  "charging”  jog  see  text.) 


Only  seven  ol  the  eight  samples  were  deemed  usable.  Two  hundred  and  thirty  two 
sets  ol  three  parallel  lines  (0. 1 2 cm  long)  were  etched  into  the  oxide  layer  of  the  samples, 
Figure  14  presents  a map  depicting  the  areas  of  pattern  delineation  on  one  of  the  samples. 
Photomicrographs  of  typical  line  sets  are  shown  in  figures  I 5 through  1 7. 


Figure  14.  Map  for  sample  I from  group  II, 
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l igurc  15.  Rellectum  pliotmiiicioiii .ipli  ( 5()()\ ) ol  sample  I ol  fiioup  1 . (;i > .loji  in  lino  set  1 probably  caused 
by  sample  charging,  (b)  I oss  ol  patient  definition  caused  by  excessive  undeicutting  during  etching  (line  set 
2).  (c)  l ine  broadening  called  b\  uiuleicul ting  (line  set  4).  (d)  better  contiol  ol  etching  (line  se  )• 


F 


f 

| 


l iga.e  16.  Reflection  photomicrograph  (500X)  ol  sample  2 of  group  II.  Lines  pass  through  many  surface 
irregularities  (discussed  in  text),  (’holograph  is  typical  of  all  line  sets  on  this  sample. 
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f igure  1 7.  Typical  1000X  photomicrographs  of  etched  line  sets  on  sample  1 of  group  II 
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Figures  1 3 and  1 5(a)  reveal  an  early  problem  that  does  not  significantly  affect  light 
scattering  in  the  guides.  The  jog  in  the  line  set  may  have  been  caused  by  stage  translation 
during  exposure,  or  more  likely  by  specimen  charging,  which  momentarily  deflected  the 
electron  beam.  The  latter  cause  can  be  eliminated  by  coating  the  resist  with  a thin  electri- 
cally conductive  film. 

The  width  of  the  354  lines  in  the  59-line  sets  on  sample  1 of  group  !1  is  1.3  +0.1  jnn; 
the  length  of  the  lines  extends  the  full  crystal  width,  which  is  1 . 1 cm.  The  separations  of  the 
lines  arc  constant  within  each  line  set,  but  vary  among  the  sets  from  2.5  to  4.0  pm  in  controlled 
steps.  An  effort  was  made  to  align  the  line  sets  for  exposure  parallel  to  most  of  the  crystal 
defects;  however,  in  some  cases  the  lines  unavoidably  intersect  crystal  defects.  Table  ? lists 
the  values  of  the  A and  B dimensions  for  the  59  sets  of  six  lines.  Dimension  A was  fixed  at 
1 .3  pni  and  is  seen  to  have  come  out  to  be  1 .3  i0. 1 pm  after  the  pattern  definition  and 
etching  processes  were  completed.  Dimension  B was  varied  from  2.5  to  4.0  pm  in  approxi 
rnately  0.3-pm  steps,  with  several  groups  of  sets  made  at  the  narrowest  dimension  that  came 
out  to  have  2,5-,  2.6-,  and  2.7-pm  separations.  The  other  steps  turned  out  to  be  2.9,  3.3,  3.6, 
and  4.0  pm. 


TABLE  2. 

LINE  DIMENSIONS  FOR  SAMPLE  1 OF  GROUP  II. 


29 


V 


LASER  DEVELOPMENT 


INTRODUCTION 

Work  lias  continued  at  RCA  Laboratories  and  California  Institute  of  Technology  on 
developing  light  sources  for  IOC  and  Fiber-Optic  applications.  The  work  tails  into  two  cate- 
gories: ( 1 ) continuous  room-temperature  injection  lasers  for  single-mode  or  multimode 
libers  and  for  IOC's,  and  ( 2. ) novel  injection  laser  sources  employing  distributed  feedback 

RCA  Laboratories  has  addressed  the  first  area,  while  the  distributed-feedback  (DFB) 
injection  laser  work  has  been  undertaken  by  California  Institute  ot  Technology. 

CW  LASER  DIODE  RESEARCH 

Excellent  control  of  the  various  parameters  involved  in  the  fabrication  of  cw  lasers 
has  been  achieved.  Some  half-dozen  wafers  have  been  grown,  all  yielding  cw  material.  The 
effect  of  diode  parameter  variations  on  device  performance  has  been  determined  in  a number 
of  cases,  and  acceptable  limits  have  been  established.  Experiments  and  studies  were  made, 
using  diodes  developed  in  this  program,  concerning  the  coupling  and  transmission  problems 
arising  in  fiber  communication  systems. 


MATERIALS  TECHNOLOGY 

The  Dll  structure  with  A1  in  the  active  region  is  now  being  used  consistently.  The 
temperature  dependence  of  the  threshold  current  density  is  relatively  steep  in  these  struc- 
tures, as  shown  in  figure  18.  This  is  partly  due  to  the  us^  of  an  undoped  active  layer  but  may 
also  be  due  to  the  A1  content  of  the  active  layer.  Thus,  experiments  are  needed  in  which  the 
threshold  dependence  is  evaluated  for  doped  active  layers,  as  it  is  possible  that  lower  threshold 
currents  may  be  achieved  thereby. 


DEVICE  DESIGN 

Careful  studies  have  been  made  of  the  temperature  increase  due  to  cw  operation.  It 
has  been  concluded  that  the  first  GaAlAs  layer  bounding  the  active  layer  exerts  a crucial 
role  in  determining  the  junction  temperature  rise  because  of  two  effects:  lirst.  the  thermal 
conductivity  of  GaAlAs  compared  with  that  of  GaAs  is  reduced,  and,  second,  the  stripe 
geometry  does  not  affect  regions  adjacent  to  or  coinciding  with  the  heat  source.  This  is 
illustrated  in  a calculation  summarized  in  figure  19.  In  other  words,  the  stripe  geometry 
reduces  the  effect  of  thermal  or  electrical  barrieis  some  distance  from  the  active  layer  (for 
example,  the  region  near  the  contacts',  but  it  does  not  alter  the  thermal  resistance  of  the 
GaAlAs  layer.  Thus,  a careful  control  of  that  layer's  thickness  becomes  important;  however, 
one  cannot  reduce  the  thickness  too  much,  because  the  optical  field  spreads  through  this 
layer  until  it  reaches  the  highly  absorbing  GaAs  regions.  An  excessively  thin  GaAlAs  layer 
would  lead  to  high  absorption  losses  and  increased  thresholds. 

The  above  problems  are  intimately  connected  with  the  question  ot  beam  width. 
Allowing  more  of  the  energy  to  leak  out  of  the  active  region  by  reducing  the  heterojunction 
barrier  height  leads  to  both  reduced  beam  widths  and  higher  thresholds.  The  best  units  have 
threshold  current  densities  on  the  order  of  1 100  - 1300  A/cm-  and  55°  beam  widths  (full  ^ 
width  at  1 /2  power  in  a direction  normal  to  the  junction).  Reducing  the  beam  width  to  35 
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STRIPE 


FOR  W 100  # # 

AT  1.6°  (A/GaAs)  + 0.4  (GaAsI  + 4 (In) 
♦ 22#(Cu)  8°  ABOVE  TQ 


AT  = + HS(W>] 

FOR  V\i  20jjm  0 

AT  - 1 .5  (A  IGaAs)  + 0.32  (GaAs)  + 2.5  (In) 
+0.58°  (Cu)  = 4.9°  ABOVE  Tq 


Figure  I V.  Calculation  of  icmpcraiuro  drops  in  the  one-dimensional  (broad  arcc) 
and  I lie  two-dimensional  (stripe)  geometries. 


tas  increased  Hk  threshold  by  a factor  of  about  2.  which  is  more  than  can  be  tolerated.  In 
view  of  the  importance  of  the  bcamwidth  in  potential  diode  applications,  a considerable 
effort  will  be  devoted  to  the  study  of  this  problem. 


1)101)1  ASSEMBLY 

Attention  has  been  focused  recently  on  a technological  problem  in  the  utilization  of 
cw  lasers,  a problem  that  arises  because  of  two  conflicting  requirements.  On  the  one  hand 
the  laser  is  fragile  and  sensitive  to  ambients,  so  that  one  would  want  to  enclose  it  in  a sealed 
structure  with  .suitable  windows  for  the  emission.  On  the  other  hand,  users  of  these  units,  a 
the  present  state  of  development,  desire  to  experiment  with  the  coupling,  and  therefore 
prefer  to  have  the  facet  accessible.  A new  design  has  been  developed  which,  to  some  extent, 
satisfies  these  requirements.  The  laser  diode  is  mounted  in  a deep  groove  milled  in  a copper 
stud  with  parts  of  the  stud  extending  and  protecting  the  diode,  figure  20  shows  this  package 
mounted  in  a liquid-cooled  heat  sink.  Ambient  protection  is  at  present  provided  by  external 
air  control,  but  may  ultimately  be  provided  by  evaporating  a suitable  oxide  film  over  le 
facet  of  the  whole  diode.  The  water  cooling  will  be  replaced  by  thermoelectric  cooling  in 
practical  mounts.  Another  method  of  laser  mounting  uses  the  “L-mount  assembly  shown 

in  figure  21 . 

1)101)1  l-'IBI  R COUPLING* 

Since  fiber  optical  communications  links  are  often  limited  in  their  performance  by 
the  amount  of  light  reaching  the  detector,  it  is  important  to  couple  into  the  fiber  as  much  of 
the  light  emitted  by  the  source  as  possible.  Several  cases  are  of  interest:  LLD  and  injection 
laser  sources,  single  fibers  and  bundles  of  fibers,  and  large-aperture  and  small-aperture  types 
of  fibers  These  variables  can  occur  in  various  different  combinations  in  different  propose 
systems.  To  provide  insight  into  the  problem  of  how  to  optimize  coupling,  measurements 

*Ttus  work  ai  RCA  was  sponsored  by  the  Ot t ice  of  Naval  Research. 
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were  made  on  the  beam  patterns  of  “normal”  rectangular  junction  plane  LLDs,  on  “stripe- 
contact”  LH)s.  and  on  a ew  injection  laser.  In  all  cases,  only  an  edge-emission  geometry 
was  considered,  as  it  is  believed  that  this  geometry  will  prove  most  practical  in  actual 
systems. 

The  beam  patterns  of  the  LL  Ds  show  slight  differences  with  polarization,  but  these 
are  small  enough  that  the  incoherent  emission  can  be  considered  unpolarized.  This  is  not 
true  of  the  laser  emission,  which  was  about  90/?  polarized  when  operated  10  - 20 7,  above 
threshold. 

In  general,  the  emission  pattern  in  the  plane  of  the  p-n  junction  differs  from  that  in 
the  plane  perpendicular  to  the  junction,  anu  the  laser  pattern  differs  significantly  from  the 
LI  D patterns.  1 or  all  the  incoherent  LLDs,  the  emission  pattern  in  the  plane  of  the  p-n 
junction  approximates  the  cos  0 pattern  associated  with  a Lambertian  source.  However,  in  the 
plane  perpendicular  to  the  junction,  the  beam  patterns  differ  significantly  from  a cos  0 
dependence  and  from  each  other,  depending  on  the  details  of  the  epitaxial  layers  and  their 
compositions  and  thicknesses  (see  figures  22  and  23). 

It  is  a good  approximation  to  assume  that  the  relative  intensity  at  any  polar  angle 
(0.  0)  from  the  axis  of  the  device  is  given  by  the  product  of  the  relative  intensities  in  the 
corresponding  angles  in  these  two  planes.  With  this  assumption,  the  relative  intensity  distri- 
bution, averaged  over  the  azimuthal  angle  0,  can  be  computed  and  integrated  over  the  solid 
angle  about  the  system  axis.  This  yields  the  cumulative  distributions  shown  in  figure  24, 
along  with  the  corresponding  Lambertian  source  distribution.  It  is  seen  that  the  beam  patterns 
from  the  LEDs  are  somewhat  more  forward-directed  than  they  would  be  from  a Lambertian 
source,  thus  showing  the  (weak)  optical  guiding  of  the  epitaxial  layers.  For  low  cone  angles, 
as  used  in  coupling  an  LED  to  a small-aperture  fiber  by  “butting”  the  filler  end  against  the 
emitting  facet  of  the  LED,  this  “beam  narrowing”  can  result  in  twice  as  much  light  coupled 
into  the  liber  as  would  be  the  case  from  a Lambertian  source.  Less-directional  LEDs  (for 
example,  curve  (b),  in  figure  24)  give  considerably  less  improvement.  In  any  case,  only  about 
47/  (—14  dB ) of  the  light  emitted  from  the  (one)  facet  can  be  coupled  into  an  N.A.  = 0.14 
fiber  in  this  way. 

Various  techniques  have  been  proposed  to  increase  the  coupling,  by  utilizing  the 
light  emitted  into  a much  larger  cone  angle;  this  can  be  done  with  either  lenses  or  reflectors. 
Well  known  principles  of  optics  prevent  this  from  being  done  without  enlarging  the  effective 
source  size,  however,  The  coupling  of  an  edge  emitter,  whose  junction  width  is  (or  can  be 
made)  greater  than  the  diameter  of  a fiber,  cannot  be  usefully  increased  in  this  one  dimension 
by  any  such  optical  elements.  Perpendicular  to  the  plane  of  the  junction,  where  the  thick- 
ness of  the  emitting  region  is  only  a few  micrometers,  the  light  emitted  into  larger  beam 
angles  than  those  corresponding  to  the  N.A.  of  the  fiber  can  be  collected  with  lenses,  etc, 
and  brought  within  the  acceptance  cone  of  the  fiber. 

For  best  effect,  lenses  used  for  this  purpose  must  be  very  small.  If  cylindrical  lenses 
are  used,  as  the  above  considerations  suggest,  the  lack  of  focusing  in  the  plane  of  the  junction 
requires  that  the  fiber  end  be  as  close  as  possible  to  the  LED.  1 his  requires  small  optical 
elements.  If,  instead,  spherical  lenses  are  used,  the  strong  aberrations  associated  with  this 
geometry  prevent  improved  coupling  unless  the  lenses  are  again  comparable  in  size  to  the 
fiber  diameter.  Using  such  small  elements,  we  have  achieved  a factor-ol-2  improvement  in 
light  coupling  over  simple  butting  in  certain  LEDs.  However,  for  coupling  into  single  fibers, 
alignment  tolerances  are  somewhat  critical. 
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Figure  24.  Fraction  of  emitted  beam  witliin  a cone  of  angle  0 (centered  on  the  device  axis)  for  two  LEDs 
and  an  injection  laser  operating  10  - 20%  above  threshold.  LED  curve  (a)  corresponds  to  the  LED  indicated 
by  the  (•)  points  of  figures  22  and  23,  curve  (b)  to  the  (o)  points.  The  corresponding  Lambertian  cumulative 
distribution  is  also  shown. 

The  situation  is  considerably  improved  if  one  is  coupling  into  a bundle  of  fibers. 

Now,  despite  losses  associated  with  the  “packing  fraction”  of  the  bundle,  source  magnifica- 
tion in  both  directions  can  be  tolerated  and  significant  improvements  can  be  achieved  by  use 
of  focusing  elements.  However,  the  improvement  to  be  expected  even  here  is  reduced  by 
reflection  and/or  transmission  losses. 

Present  cw  room-temperature  injection  lasers  are  typically  multimode  devices  when 
operated  10  — 20%  above  threshold.  It  is  thus  more  difficult  to  couple  the  output  beam  into 
a single  fiber  than  would  otherwise  be  the  case.  A typical  laser  has  the  beam  pattern  shown 
in  figure  25.  The  calculation  outlined  above  for  LEDs,  when  applied  to  this  laser,  gives  the 
results  shown  in  figure  24.  It  is  seen  that  the  much  more  directive  beam  pattern  of  the  laser 
permits  approximately  1 5%  of  the  output  to  be  coupled  into  the  fiber  with  simple  butting. 
Simple  (small)  cylindrical  lenses  permit  this  to  be  significantly  increased,  with  increased  con- 
struction and  alignment  complexity. 

The  beam  pattern  of  an  injection  laser  varies  with  excitation  level,  from  the  broad 
LED-type  pattern  below  threshold,  through  the  “beam-narrowing”  region  of  increasing 
stimulated  emission  as  threshold  is  reached,  to  the  quite  sharp  pattern  slightly  above  threshold 
when  only  one  inode  (or  a few)  is  oscillating.  At  the  still  higher  excitation  levels  needed  to 
provide  optimal  power  output,  the  number  of  laser  modes  oscillating  increases  sharply,  and 
the  beam  pattern  again  broadens.  Coupling  efficiency  is,  therefore,  a somewhat  ill-defined 
term  for  these  devices. 

This  characteristic  is  illustrated  in  figure  26,  where  the  coupling  loss,  measured  as  the 
ratio  of  power  coupled  through  1000  feet  of  N.A.  = 0.14,  low-loss,  step-index  fiber  to  the 
power  measured  into  a large-area  detector  butted  against  the  emitting  facet  of  the  laser,  is 
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plotted  against  laser  drive  current,  as  is  the  power  through  the  fiber.  The  figure  shows  a 
gradual  increase  in  coupling  efficiency  as  threshold  is  approached,  with  a sharp  rise  just  above 
threshold.  At  still  higher  drive  levels,  the  beam  pattern  broadens  and  the  coupling  efficiency 
drops  somewhat. 

LIFE  TESTS 

Excellent  progress  is  being  made  in  eliminating  the  causes  of  rapid  diode  failure.  To 
date  the  longest  operating  time  in  the  testing  laboratory  is  2000  hours.  Figure  27  shows  the 
power  output  vs  current  curve  of  a diode  randomly  selected  from  a wafer  at  the  start  of  the 
test  and  alter  2000  hours  of  operation.  The  threshold  current  increase  is  about  10'#.  There- 
fore, the  power  output  change  with  time  will  depend  on  the  operating  current  above  threshold 
Unfortunately,  a power  failure  in  the  testing  laboratory,  with  a resultant  “surge,”  destroyed 
all  the  diodes  life-tested,  and  no  additional  time  substantially  in  excess  of  2000  hours  could 
be  accumulated.  New  facilities  have  been  constructed  with  a special  protection  circuit  which 
will,  in  the  future,  shut  the  diodes  off  permanently  in  case  of  power  failure.  More  diodes 
are  being  tested  and  data  will  be  accumulated. 

The  degradation  mechanism  operating  in  these  units  is  predominantly  facet  damage. 
Presumably,  portions  of  the  material  lose  their  volatile  component  (As),  leaving  a residue  of 
Al  and  (la.  Early  results  indicate  that  this  problem  can  be  drastically  reduced  through  the 
deposition  of  a protecting  AI2O3  film  on  the  facet.  Experiments  are  continuing  toward 
improved  deposition  techniques  as  the  device  yield  at  present  is  still  not  satisfactory. 


Figure  27.  Power  output  as  a function  of  de  current  of  a diode  operating  at  room  temperature  before  and  after 
2000  hours  of  continuous  operation  in  dry  air  ambient. 
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DISTRIBUTED-FEEDBACK  L\SERS 

Prompted  In  recent  experimental  * and  theoretical--  work,  there  has  Been  much 
interest  in  applying  the  concept  of  distributed  leed back  to  (iaAs  lasers  where  control  ol  the 
emission  spectrum  is  difficult.  Shank  et  al--*  have  demonstrated  that  one  can  reduce  the 
threshold  ol  an  optically  pumped  distributed-leedhack  <1)1  B)  laser  by  using  a double  hetero- 
junction (DID  rather  than  a single  heterojunction  (SID.  This  demonstration  points  the  way 
to  the  ultimate  realization  of  a 1)1  B injection  laser.  Moving  in  this  direction.  Scifres  et  al- 
constructcd  a 1)1  B single  heterojunction  (SID  (iaAs  diode  laser  which  was  electrically  pumped 
and  had  threshold  current  densities  comparable  to  those  of  normal  SH  diodes.  The  following 
section  is  a summary  ot  the  work  done  by  A.  ^ .'iriv  and  his  coworkers  at  ( alilornia  Institute  ol 
Technology  and  \1  Nakamura--''  of  the  Hitachi  Central  Research  Laboratories  in  this  area  during 
this  quarter. 

THRESHOLD  SHIFTS  IN  Gaj  Al\  AS  DOUBLE  HETEROJ UNCTION  LASERS 
WITH  A CORRUGATED  INTERFACE26 

Consider  the  double  heterojunction  structure  in  figure  28.  The  oscillating  optical 
fields  associated  with  corrugated  structure  obey  the  following  coupled  relations: 

— = iK^  exp  |-2i6,  /.}  A+  aA* 
dz 

and 

= i Ky  exp  1 2i5^  /.|  A + aA+. 
d/ 


where 

A+.  A"  = right- and  left-going  normal  mode  amplitudes,  respectively 

A = corrugation  period 
t!  ^ re  licet  ion  order 

= propagation  constant  ol  the  waveguide  layer  in 
2 a = intensity  gain  constant  of  the  modes 

= coupling  constant  (rellection  produced  by  the  corrugations) 

Kp'or  the  m1'1  order  longitudinal  mode  of  rectangular  corrugations  is  given  by 

-i  v ■>  2.  1 
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In  the  limit  X/a  >>  I that  is.  the  corrugation  height  being  much  smaller  than  the 
wavelength  the  above  reduces  to 


(-1  l^Xa 
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lor  the  lowest  order  whieh  is  independent  of  n, . Thus,  a corrugation  height  of  a is  equa  y 
effective  at  either  an  air-T.aAs  or  C.aAs  -C,a,.x  AlxAs  interlace.  Fins  raises  the  poss.l  y 
of  using  internal  corrugations  to  achieve  feedback  where  one  might  expect  that  tl« : rcllcct 
coupling  would  be  severely  reduced  due  to  the  small  index  difference  available,  gurcs- 
and  30  show  numerical  solutions  for  the  threshold  gam  ol  a corrugated  inte'Hace  hctc 
junction  laser  calculated  from  the  reflection  coupling  constant  kfc  with  C ~ 1 • A,  ° 
plotted  in  each  figure  are  threshold  gam  curves  lor  a C.aAs  laser  with  cleaved  y 
reflection  end  mirrors.  Note  that  for  L < 1 nun  and  t < 0.75  gm.  the  com.ga  ed  mterf 
lasing  threshold  is  much  lower  than  that  of  the  Fabry-Perot  cavity.  One  can  also  sec  t 
the  transverse  mode  number  is  dependent  on  both  the  threshold  and  the  indices  n,  and  n2. 
Thus  control  of  the  coupling  constant  will  also  control  the  transverse  mode. 
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FABRICATION  AND  ANALYSIS  OF  INPUT/OUTPUT 
WAVEGUIDE  COUPLING  STRUCTURES 
(UNIVERSITY  OF  WASHINGTON) 


FIBER  COUPLING  TRANSITION  ELEMENTS 

The  general  problem  of  coupling  from  integrated  optical  circuits  to  fiber  transmission 
lines  is  dominated  by  practical  considerations  of  the  small  dimensions  and  installation  difficulties 
involved  in  implementing  any  particular  scheme.  A useful  method  for  coping  with  these  prob- 
lems involves  employment  of  transition  fibers  in  the  manner  illustrated  in  figure  31.  The  transi- 
tion fiber  is  coupled  on  one  end  to  the  output  ot  an  integrated  optical  waveguide  and  on  the 
other  end  to  a fiber  transmission  line.  The  former  coupling  is  effected  under  controlled  condi- 
tions with  accurate  high-resolution  equipment  that  represents  part  of  the  overall  integrated 
optical  circuit  fabrication  process.  The  latter  coupling  must  be  effected  in  the  field  with  jigs 
of  the  kind  that  will  be  employed  for  general  fiber  splicing.  The  transition  fiber  has  a length 
of  perhaps  1/3  m 

A practical  method  for  implementing  the  transition  fiber  involves  use  of  an  externally 
mounted  fiber.  This  is  a fiber  mounted  on  the  surface  of  a larger  support  fiber  in  the  manner 
shown  in  figure  32.  The  externally  mounted  fiber  is  distributive^  coupled  to  the  waveguide 
at  one  end  and  butt-joined  to  the  fiber  transmission  line  at  the  other.  Distributed  coupling 
is  most  readily  effected  when  the  externally  mounted  fiber  has  a rectangular  cross  section. 

The  short  thickness  dimension  is  on  the  order  of  1 micrometer  to  permit  only  single-mode 
operation  in  this  direction.  The  width  lies  in  the  range  of  5 - 30  micrometers  in  accordance 
with  the  taper  that  can  be  provided  in  the  integrated  optical  waveguide.  The  support  fiber 
lias  a width  of  50  - 100  pin.  This  dimension  is  established  by  considerations  ot  the  con- 
venience with  which  the  waveguides  can  be  separated  on  the  integrated  optical  substrate. 


figure  31.  lOC-to-fiber  coupling. 


Distributed  coupling  can  be  provided  with  the  aid  of  tapered  films  or  variable-period 
periodic  structures  in  the  manner  illustrated  in  figure  33.  The  effectiveness  of  these  structures 
has  been  established  on  an  analytic  basis,  and  further  research  will  be  directed  toward  imple- 
menting these  structures  and  evaluating  potential  sources  of  practical  difficulty.  The  butt- 
joined  end  of  the  taper  will  require  a circm.ir  output  end  for  the  transition  fiber.  This  may 
be  accomplished  by  altering  the  cross  section  of  the  externally  mounted  fiber  before  or  during 
the  fabrication  process  as  indicated  in  figure  35  or  by  providing  a second  coupling  region. 
Figure  34  illustrates  how  the  second  coupling  step  may  be  accomplished  with  a second  extern- 
ally mounted  fiber  positioned  either  on  a separate  support  fibir  or  on  the  same  one.  1 ollowing 
fabrication  of  the  coupling  regions,  the  externally  mounted  fibers  will  be  coated  for  protection 
against  dust  and  other  sources  of  optical  transmission  loss. 

Two  topics  relevant  to  implementation  of  the  coupling  scheme  outlined  above  are 
discussed  in  this  section.  The  first  concerns  the  solution  of  coupled-mode  equations  when  the 
propagation  constant  in  two  coupled  waveguides  varies  linearly  with  opposite  slope.  The  maj  'r 
result  is  a relatively  simple  asymptotic  form  for  the  power  coupled  from  one  guide  to  another 
This  expression  shows  that  virtually  complete  coupling  occurs  so  long  as  the  coupling  coeffi- 
cient exceeds  some  minimum  value  and  so  long  as  the  coupling  region  is  of  sufficient  length. 
Exact  knowledge  of  the  average  mode  propagation  constants  is  not  required.  The  second 
topic  is  a review  of  the  experimental  results  achieved  so  far.  Externally  mounted  fibers  have 
been  fabricated  and  preliminary  coupling  measurements  have  been  performed. 


DESIGN  OF  DISTRIBUTED  COUPLING  STRUCTURES 


Apart  from  the  practical  problems  of  fabricating  coupling  structures  between  optical 
guides,  there  are  conceptual  design  questions  which  deserve  study.  Even  if  the  coupling  and 
propagation  parameters  can  be  controlled  approximately  in  the  fabrication  process,  there 
exists  the  question  of  what  the  parameter  values  should  be  in  order  to  achieve  efficient  coup- 
lers with  predictable  performance.  Two  aspects  of  this  design  question  have  been  studied; 
the  first  concerns  minimizing  the  effect  of  uncertainties  and  variations  in  parameter  values  on 
coupler  performance,  and  the  second  concerns  optimizing  the  parameters  when  their  values 
can  be  accurately  controlled.  Achieving  an  effective  coupler  in  a minimum  length  is  an  example 
of  a design  objective. 

Under  assumptions  which  are  widely  applicable,  the  coupling  of  power  between  guided 
modes  can  be  described  in  terms  of  the  so-called  coupled-mode  equations,  and  this  formalism 
has  been  the  subject  of  the  present  study.  Particular  attention  has  been  paid  to  obtaining 
coupler  designs  which  are  insensitive  to  the  parameters  of  the  fabricated  device. 


SOLUTIONS  OF  THE  COUPLED-MODF  EQUATIONS 

The  particular  model  considered  is  described  by  the  coupled-mode  equations  in  vector 
form  for  propagation  in  the  z variable: 


da(z)  _ . 
dz 
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where  a](z)  and  aifz)  are  the  individual  mode  amplitudes,  and  k | fz).  lofz),  cj c'2J  are  the 
respective  propagation  constants  and  coupling  coefficients.  A constant  complex  coupling 
coefficient  is  assumed,  c 1 2=c2 1 =c.  It  is  also  necessary  to  assume  that  |c|  « kj,  k^.  The 
k(z)  are  taken  to  vary  linearly  at  rates  set  by  a parameter  V>  0 introduced  through  the 
relations 

k|(z)  = kQ(  I+7Z),  and  k^fz)  = k0(  I-7Z)  . ]0) 

Exact  phase  match  occurs  at  z = 0.  where  both  propagation  constants  equal  k0.  Approximate 
phase  match  and  significant  power  transfer  occur  over  a distance  inversely  proportional  to 
y and  centered  at  z = 0.  For  sufficiently  large  z.  say  for  0.2  < |yz|,  the  phase  mismatch 
prevents  power  exchange  between  modes.  In  physical  couplers  the  point  z = 0 may  not  be 
known,  but  this  introduces  effects  that  are  not  central  to  the  present  analysis. 

The  differential  equations  are  further  reduced  by  treating  the  slowly  varying  mode 
amplitudes  A(z.)  by  writing  a(z)  as 

a(z)  = exp(ikt)z)  • M /)  . (11) 
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The  independent  variable  z is  reduced  to  dimensionless  form  by  defining 


x = s/27k0  • z . 

To  distinguish  the  mode  amplitudes  as  they  are  dependent  on  z or  x,  we  use  the 
notation  B(x)  with  the  equivalences 

A(z)  h B(s/27k0  L)  = B(x)  . 


(12) 


(13) 


With  these  changes  of  variables  the  coupled  first-order  equations  become 
x/2  c/\f2ykQ 


d B(x) _ . 
dx 


B(x) 


(14) 


c*lJ2ik~0  - x/2 

Further  differentiation  yields  the  separated  second-order  equations: 
d"Bi(x) 

+ [x2/4  + k -i/2]  B,(x)  = 0 


and 


dx* 


d“Bi(x) 

+[x2/4  + ^ + i/2]  B2(x)  = 0 

dx~ 


(15) 


where  v = |c|2/2',  kQ  is  an  important  parameter  which  deter  lines  the  *'orm  of  the  solutions. 
The  second-order  equations  are  in  a standardized  form  of  the  equation  for  the  parabolic 
cylinder  functions  22  which  is 


d2w(x) 

dx2 


+ 1 , r - a]  w(x)  = 0 . 


(16) 


For  the  present  problem  the  upper  sign  applies,  and  the  parameter,  a,  takes  on  a different 
value  for  each  of  the  two  modes;  denoting  these  as  aj  and  a'2,  respectively  (not  the  mode 
amplitudes),  we  have 


a',  = - v + \/2  a2  = -v  - 'ill  . (17) 

The  mathematical  problem  now  is  to  find  solutions  for  B](x)  and  B2(x)  in  terms  of  the 
parabolic  cylinder  functions  for  parameter  values  a'j  and  ah.  These  solutions  must 
together  satisfy  the  coupled  first-order  equations. 

The  power  in  each  mode  is  given  by  the  squared  magnitude  of  any  of  the  mode 
amplitudes.  For  real  k\  and  k2  and  for  the  stated  forms  of  the  coupling  coefficients  the 
total  power  will  be  constant.  That  is, 

Pj(z)  = lajl2  — |Ajl2  = |Bj|2  i=l,2  (18) 
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and 


P|(/)  + P->(z)  = constant 


(19) 


For  an  initial  condition  the  usual  assumption  is  made  that  all  the  power  is  in  one  mode. 
Specifically,  let  P|(z-*-°°)  = 1 and  P2(z.->-«>)  = 0.  The  corresponding  phase  of  the  mode 
amplitude  may  be  arbitrary. 

Of  the  parabolic  cylinder  functions  standardized  by  Miller, 27  the  ones  appropriate 
to  this  problem  are  the  two  independent  solutions  he  denotes  by  E(a,x)  and  E*(a.x).27 
The  properties  listed  in  reference  27  for  these  functions  assume  that  the  parameter  a is  purely 
real.  For  the  present  problem  both  values,  aj  and  a2,  are  complex  and  some  of  the  listed 
expressions  cannot  be  used  without  being  generalized  to  the  case  of  complex  a.  This  gen- 
eralization can  be  effected,  however,  ai:d  for  the  specific  values  of  a’]  and  aS  these  solutions 
may  be  denoted  by 

Fj(x)  = H(aj,x)  and  E*(x)  = F*(aj,x)  i = 1 , 2 (20) 

With  these  functions  the  solutions  for  the  mode  amplitudes  have  been  found  to  be 
A | (z)  = Bj(x)  = b*H|(-x) 

= b’(o)s  X * Li  | ( x ) + i sin  A hj(x)l 

and  (21) 

A2(z)  = B2(x)=b*E2'-x) 

= b*|cos  X E2(x)  - i sin  X E*(x)]  , 

where  the  coefficient  b = [i>'/2  exp(i^)/2] 1 - and  ip  is  the  phase  of  the  coupling  coefficient, 

— ph  c.  The  term  X is  used  for  convenience  in  writing  some  coefficients  depending  on  v 
and  may  be  considered  to  be  defined  by  the  relations 

cos  A = exp(-7n>)  and  sin  X = ( 1 - exp(-2in>)) ' 2 . (22) 

Also  recall  that  v=  |c|“/27V0  and  x =s]2yk0  z.  The  second  form  of  these  solutions 
(involving  X in  (21 ) results  from  the  fact  that  E(a,-x)  may  be  written  as  a linear  combina- 
tion of  E(a,x)  and  E*(a,x).  These  second  forms  are  preferred  for  positive  x because  the 
available  asymptotic  expressions  for  E(a,x)  are  valid  only  when  x is  positive.  It  is  these 
asymptotic  expressions  which  make  it  possible  to  obtain  the  simple  expressions  for  coupler 
performance. 


CHARACTERISTICS  OF  THE  SOLUTION 

Writing  the  power  as  the  mean  squared  magnitude  of  the  given  solutions  of  the 
coupled-mode  equations  and  using  the  asymptotic  expressions  for  the  E(a,x)  functions, 
we  obtain  the  following  results  for  large  x (neglecting  terms  of  higher  order  than  1/x-). 

Pj(x)  = |B|  (x)|“  'V  | - !>/x“  ->  1 o«  X“,  x < 0 


and 


(23) 


For  the  case  of  positive  x we  find: 

Pj(x)  ^ eos~X  + (2v  2/x)sinXcosXsin[ x“/2  + 2t^fenx  + <p  + 7r/4) 


and 


+ (t>/x~)(sin”X  - cos~X) 


Pifx)  'v  sin “X  ~(2v  2/x)sinXcosXsin[  x“/2  + 2t*Cnx  + 0 + 7r/4  j 


- (»*/x“)(sin“X  - cos“X)  (24) 

The  condition  for  validity  of  these  forms  is  that  s/l+4v-  «x2.  the  coefficients  in  X have 
been  given  above,  and  0 is  a phase  angle  varying  only  slowly  with  v. 

In  (24)  the  leading  term  is  a constant  which  is  the  power  remaining  in  each  mode  as 
x -*  00.  That  is: 

P|(+oo)  = exp(-2w),  |>2(+oo)  = 1 _eXp(-2w),  (25) 

and  p2(+°o)  is,  ot  course,  just  the  coupler  efficiency  provided  the  coupling  region  is  long 
enough.  The  record  term  represents  an  oscillatory  exchange  of  power  between  the  modes 
which  decreases  in  magnitude  as  l/x.  The  frequency  of  oscillation  increases  without  limit  as 
x/2  (representing  the  increasing  phase  mismatch):  and  the  magnitude  of  the  oscillation,  at  any 
fixed  x,  is  proportional  to  2 /<  2 exp(-Trt>)  [ I - exp(-2rrt*)]  which  decreases  exponentially  for 
large  v.  These  general  features  ol  the  power,  in  mode  2 say,  are  shown  in  figure  35,  which  is 
drawn  for  v ^ 0.25.  That  this  concept  could  lead  to  high  efficiency  for  coupling  to  the 
second  mode  was  well  known  from  the  microwave  ease,  but  the  simple  expression  for 
predicting  this  efficiency  from  the  parameters  is  new. 

A program  for  the  direct  numerical  integration  of  the  coupled-mode  equations  has 
been  written  and  used  on  several  cases  representing  this  coupler.  The  results  are  shown  in 
figures  36  38,  and  fully  confirm  the  properties  found  analytically.  In  particular,  the 
power  approaches  the  constant  value  predicted  in  each  case. 

The  expression  for  the  asymptotic  efficiency  of  coupling  from  the  first  to  the  second 
mode  is  remarkably  simple,  and  has  a simple  dependence  on  the  parameters.  It  is 

FIT  - I - exp(-27rr')  (26) 

with 

v — |c | ~/2"yk0  = - lc/kQ|“  • (kQ/7)  . (27) 

The  efficiency  depends  on  c and  7 only  through  the  value  of  v.  Couplers  with  differing 
values  of  c and  7 can  have  the  same  efficiency  if  they  have  the  same  |c|  ~/y  ratio.  The 
second  form  for  v above  shows  the  convenience  of  expressing  both  c and  7 in  units  of 
k0.  Represent:* live  values  of  v are  0.1  for  an  efficiency  of  50%  and  0.50  for  96%.  For 
even  greater  values  oft*  the  efficiency  becomes  arbitrarily  close  to  unity,  and  it  becomes 
increasingly  insensitive  to  uncertainties  in  the  values  of  e and  7. 

It  is  convenient  to  show  these  results  graphically.  This  is  done  in  figure  39,  where 
the  values  of  c and  7 necessary  to  obtain  given  efficiencies  are  plotted.  It  should  be  recalled 
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that  the  eon  pled -mode  theory  assumes  |e| « kQ,  say  |e|  /kG  < 10"“.  Consequently,  for 
high  efficiency,  the  range  of  useful  7 will  be  10-4  ^7/k0,  efficiency  increasing  for 
decreasing  values  of  7. 

Apart  from  the  value  of  i> , which  affects  only  the  asymptotic  efficiency,  the  value 
of  7 sets  the  scale  factor  between  the  / distance  and  the  x variable:  it  will  determine  the 
necessary  length  of  a coupling  region.  Smaller  values  of  7 (chosen  for  increasing  the  efficiency) 
will  require  larger  coupling  lengths.  The  coupler  efficiency  will  be  independent  of  the  exact 
length  only  if  it  is  sufficient  for  the  oscillatory  exchange  of  power  to  have  decreased  to  some 
“small”  magntiude.  For  example,  if  v = 0.5  and  the  corresponding  efficiency  is  96%.  then  the 
power  oscillation  between  modes  becomes  less  than  ±4%  (that  is,  1-0.96)  for  x-  550.  In 
terms  of  the  physical  distance  z,  this  condition  is  that  kGz  55(k0/7  )'  2 = 5 k0/  Id  . If  the 
particular  value  of  v = 0.5  is  achieved  by  the  choice  7 /k0  = 10'4,  |c|  /kG  = 10'-.  then  we 
need  k0z  > 500.  For  smaller  values  of  |c|  and  7,  the  coupler  will  be  correspondingly  longer. 

The  value  of  kGz  > 5k0/  |c|  may  be  compared  to  the  length  of  a phase-matched  coupler 
which  must  be  exactly  of  length  k0L  = kG7r/2  |e|  for  1007&  power  transfer.  Thus,  the  “tapered” 
coupler  will  be  three  or  more  times  the  length  of  the  phase-matched  coupler. 

EXPERIMENTATION 

Experimental  work  has  centered  on  the  development  and  testing  of  externally  mounted 
fibers  and  on  commissioning  of  the  multitarget  sputtering  facility.  Preliminary  experiments  on 
simple  waveguide  coupling  structures  have  also  been  performed. 

EXTERNALLY  MOUNTED  FIBERS 

Externally  mounted  fibers  consist  of  a rectangular  (or  other)  cross  section,  high- 
refractive-index  guiding  fiber  on  top  of  a larger,  low-index  support  fiber  (figure  40).  The 
combined  fiber  is  drawn  from  a larger  preform  of  approximately  the  same  relative  cross 
section.  In  order  for  the  fiber  drawing  process  to  be  successful,  the  two  glasses  used  for  the 
guiding  and  support  fibers  must  have  similar  mechanical  and  thermal  properties  (such  as 
softening  temperature  and  thermal  expansion  coefficient).  The  glasses  currently  used  are 
Corning  7059  (ii[)  = 1 .53)  for  the  guide  and  Corning  7740  (n^  = 1.47)  for  the  support. 

Suitably  sized  sheets  are  used  to  produce  a preform  which  is  then  drawn  down  in 
a fiber  pulling  machine  to  produce  the  final  fiber.  Typical  sizes  of  the  final  liber  are 
10  by  1 pm  for  the  guiding  fiber  and  100  by  20  jum  for  the  support. 

Although  both  guides  maintain  their  rectangular  cross  section,  some  small  changes  in 
aspect  ratio  (width  to  thickness)  have  been  observed.  In  addition,  it  has  been  found  difficult 
to  produce  guides  of  high  aspect  ratio  ( >~30: 1 ). 

Some  optical  evaluation  has  been  carried  out  on  these  fibers.  It  has  proved  possible 
to  use  prism  couplers  to  selectively  excite  the  “thickness  modes”  of  these  fiber  guides. 

Modes  of  fibers  as  narrow  as  7 /am  have  been  excited.  Because  these  guiding  fibers  have 
widths  much  greater  than  their  thickness  05: 1 ).  it  is  convenient  to  subdivide  the  modes  by 
considering  groups  of  “thickness  modes”  which  describe  the  electric  (TF)  or  magnetic  (TM) 
field  distribution  across  the  thickness  of  the  guide.  For  each  thickness  mode  there  exist  a 
number  of  transverse  modes  which  describe  the  field  distribution  across  the  width  of  the 
guide.  To  illustrate  this,  figure  41  shows  a photograph  of  the  output  of  a 7 by  2 jam  wave- 
guide coupled  through  a prism  coupler.  The  guide  was  excited  by  end  launching.  The  tour 
thickness  modes  appear  as  bands  of  mode  lines.  Each  mode  line  represents  a particular 
transverse  mode  within  that  thickness  mode. 
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COUPLING  COEFFICIENT  |c|tk 


- yi* 

(INCREASING) 


Hgurc  y>.  ( oupling  efficiency  (power  coupled  to  mode  2,  Pi) 
vs  coupling  coefficient,  c,  and  rate  of  change.  7,  of  propagation 
constants. 
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FIBER 

Figure  4U.  l:\icriKilly  mounted  fiber  consisting  of  a guiding  fiber  of 
width  w and  thickness  I on  a larger  supporl  fiber. 


TE3x  modus 


TE2x  modes 


Transmission  loss  measurements  have  been  made  for  these  waveguides  by  the 
destine  live  method  of  measuring  light  output  after  repeated  shortening  of  the  fiber.  Sub- 
strate modes  were  stripped  just  prior  to  the  fiber  end.  Losses  of  0.05  dB/cm  lor  multimode 
( 25  by  5 pin ) guide  and  0.8  dB 'em  for  guide  capable  of  carrying  only  one  thickness  mode 
have  been  measured.  J iber  drawing  has  not  as  yet  been  conducted  under  clean-room  condi- 
tions, which  can  be  expected  to  result  in  reduced  losses. 


COUPLING  I XIM  RIMI  NTS 

An  initial  investigation  ol  the  use  of  externally  mounted  fibers  as  an  interconnection 
medium  has  been  carried  out.  Preliminary  experiments  on  coupling  from  planar  guide  to 
fiber  and  from  fiber  to  fiber  have  been  performed. 


PLAN AR-GUIDl  - 1 0-1  IBI  R COUPLING.  Coupling  from  high-index  guides  on 
high-index  substrates  (such  asdittused  semiconductor  guides  or  LiNbx  Taj.x  0^  on  LiTaOO 
to  externally  mounted  libers  will  require  earelul  design  and  complicated  coupling  structures, 
since  phase  match  between  the  two  guides  cannot  be  directly  achieved.  However,  in  eases 
m which  the  planar  guide  is  ol  high  index  on  a low-index  substrate  (lower  than  the  refractive 
index  ol  the  guiding  fiber),  coupling  can  be  ellected  by  simple  tapering  of  the  planer  guide 
and  by  mounting  the  fiber  directly  over  the  taper  as  shown  in  figure  42.  Light  in  the  lowest- 
order  mode  ol  the  planar  guide  is  converted  directly  to  light  in  the  lowest-order  mode  of  the 
liber  in  this  structure,  which  represents  a type  ol  branching  waveguide  concept.  Preliminary 
experiments  indicate  that  light  is  indeed  coupled  from  the  planar  guide  to  the  fiber,  but 
accurate  ellieieney  measurements  have  not  been  performed  as  yet  because  of  dimensional 
mismatch  between  the  Gaussian  beam  propagating  in  the  planar  guide  and  the  width  of  the 
libei.  This  concept  should  be  directly  transferable  to  coupling  from  linear  guides  to  extern- 
ally mounted  fibers  of  the  same  width. 


MBI  R-TO-I  IBI  R COUPLING.  To  simulate  coupling  between  an  externally  mounted 
liber  and  a linear  guide,  experiments  were  conducted  on  coupling  between  two  externally 
mounted  libers  Irom  the  same  draw.  This  is  a special  case,  the  two  fibers  being  exactly  phase 
matched,  but  none  the  less  one  that  permits  assessment  of  the  difficulties  of  carrying  out 
the  clamping  procedure.  In  preliminary  experiments,  two  fibers  were  used  (10  by  I pm), 
with  a matching  oil  between  them  to  increase  the  coupling  coefficient.  Light  was  coupled 
into  one  liber  and  coupled  out  alter  approximately  4 cm  via  prism  couplers.  A second  fiber 
was  manipulated  into  close  proximity  of  and  parallel  to  the  first  guide.  By  changing  the 
gap  between  the  two  guides  (figure  43),  transferal  of  light  was  observed  from  the  first  fiber 
into  the  second  with  typical  efficiencies  ol  several  tens  of  percent  and  ranging  to  above  60% 

( 1007,  is  theoretically  possible).  The  experiment  as  described  is  not  sufficiently  definitive 
due  to  the  close  tolerance  or  alignment  required  (to  within  a fraction  of  a micrometer  over 
the  coupling  length)  and  not  as  yet  reliably  achieved. 
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Figure  42.  Simple  planar  guide  to  fiber  coupler. 
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Figure  43.  Schematic  of  fiberto-fiber  coupler. 


THEORETICAL  ANALYSIS  OF  COUPLING  FROM  PLANAR  TO  LINEAR 
WAVEGUIDES  VIA  HORN  STRUCTURES 

Oik*  barrier  to  efficient  coupling  front  planar  dielectric  waveguides  to  linear  waveguides 
is  the  difference  between  the  size  of  the  beam  in  a planar  waveguide  and  the  size  of  linear 
waveguides.  The  beam  width  in  planar  waveguides  is  on  the  order  of  50  - 1 00  wavelengths, 
while  the  dimensions  of  linear  waveguides  are  only  a few  wavelengths.  This  size  difference 
produces  a geometric  mismatch  whose  deleterious  coupling  effects  can  be  reduced  by  using 
a noiumiform  waveguide  or  a lens  system  to  focus  the  beam. 

There  are  many  ways  in  which  optical  coupling  between  libers  or  lasers  and  integrated 
optical  components  may  be  achieved.  In  this  section  the  use  of  a tapered  structure  similar  in 
shape  and  function  to  a microwave  horn  is  considered  as  a transition  between  a linear  wave- 
guide and  a planar  guide.  This  eollinear  coupling  technique  is  useful  under  conditions  in  which 
a beam  is  launched  in  a film  from  either  a fiber  or  a laser  and  is  then  coupled  into  a narrow 
channel  that  might  serve  as  modulator  or  switch.  The  taper  is  also  potentially  useful  tor  output 
coupling  to  fibers.  When  tapered  waveguide  section  is  used  as  a coupler,  its  efficiency  depends 
upon  the  power  that  remains  in  the  lowcst-order  mode  at  the  open  end  or  mouth  ol  the  taper. 
This  section  is  a brief  summary  of  detailed  calculations  to  be  found  in  reterence  28. 

Figure  44  shows  the  geometrical  structure  of  the  horn.  The  waveguide  considered 
has  indices:  n j = 1 .57,  m = 1 .0.  and  n4  = 1 .53.  The  thickness  in  the  y-direction  is  held 
constant  at  0,5  micrometer  and  the  thickness  tx  ol  the  linear  waveguide  is  5 micrometers. 
These  parameters  show  the  linear  waveguide  to  be  single  mode. 

The  first  shape  of  a coupling  region  to  be  investigated  is  a linear  taper.  If  the  coupling 
region  is  made  long  enough,  we  expect  to  be  able  to  couple  to  the  linear  waveguide  with  any 
degree  of  efficiency  desired;  however,  to  keep  the  structure  compact,  the  length  should  be 
kept  to  a minimum,  l or  the  linear  taper  shown  in  figure  44.  the  thickness  is  given  by 
t x = 2/tan(0),  where  0 is  the  half  angle. 

Solutions  for  the  power  contained  in  the  lowest  three  modes  in  the  reciprocal  problem 
of  coupling  from  a linear  waveguide  to  a horn  are  plotted  in  figure  45  as  a function  of  horn 
width.  These  values  are  calculated  for  an  input  power  of  I watt.  The  small  oscillations  about 
some  average  value,  observed  in  figure  45,  are  caused  by  the  coupling  of  energy  between  all 
the  modes.  This  would  disappear  if  the  horn  supported  only  two  modes. 

Beams  that  we  are  interested  in  coupling  into  linear  waveguides  contain  most  of  their 
energy  in  the  lowcst-order  mode,  at  least  when  beam  size  and  horn  mouth  are  approximately 
equal.  This  being  so.  coupling  efficiencies  can  best  be  improved  by  improving  coupling  of  the 
lowt  st-order  mode  at  the  mouth  to  the  lowcst-order  mode  in  the  waveguide.  Comparing 
efficiencies  for  coupling  to  the  lowcst-order  mode  at  difterent  angles  shows  that  care  snould 
be  taken  when  choosing  the  horn  angle;  that  is  osc'llatious  due  to  coupling  to  higher-order 
modes  may  occur  such  that  a decrease  in  efficiencies  would  result  lor  a particular  thick- 
ness although  the  horn  is  longer. 

Coupling  efficiencies  from  the  lowest-order  mode  at  the  mouth  ot  the  horn  to  a 
single-mode  linear  waveguide  can  be  improved  by  varying  the  slope  of  the  horn  wall.  Coupling 
coefficients  were  shown  to  be  proportional  to  the  slope  of  the  wall,  so  the  overall  cltect  ol 
coupling  can  be  reduced  it  dtx  <)z  is  small  when  a j /t  + as /tx  + U3/ 1^  is  large.  The  analogous 
problem  of  reflections  in  a transmission  line  has  been  investigated,  and  the  results  indicate 
this  to  be  a good  method  of  improving  efficiency.  By  allowing  the  thickness  to  vary  as 
tx  = t0  + A(  1 - c'x“/'!“  ),  the  coupling  efficiency  can  be  improved  by  10 r/r.  for  the  same- 
length  horn.  Figure  46  is  a plot  ot  power  in  the  lirst  three  modes  versus  thickness  ot  the 
horn  lor  this  thickness  variation. 
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Coupling  efficiency  for  linear  tapers  has  been  shown  to  be  approximately  50 % for 
coupling  from  a beam  whose  width  is  approximately  100  pm.  This  efficiency  is  of  course  a 
function  ot  flare  angle  and  generally  decreases  with  increasing  angle.  However,  there  is 
coupling  to  higher-order  modes,  so  there  is  some  oscillation  around  an  average  value  that 
might  tend  to  increase  the  efficiency  for  a small  increase  in  angle  at  a given  thickness. 

The  coupling  efficiency  can  be  improved  by  changmg  the  slope  of  the  waveguide 
wall  By  this  is  meant  that  the  efficiency  increases  for  a given  length  of  coupling  region, 
for  the  configuration  used  in  figure  46,  an  increase  of  approximately  10%  is  observed. 


Figure  46.  Power  in  first  three  modes  as  a function  of  horn  thickness  for  specified 
thickness  variation  (see  text). 


TRANSIENTS  IN  PERIODIC  STRUCTURES: 

COUPLED  WAVES  APPROACH 
(UNIVERSITY  OF  CALIFORNIA,  LOS  ANGELES) 

Most  of  the  work  in  the  field  of  electromagnetic  wave  propagation  in  periodic  struc- 
tures v.as  directed  toward  continuous  monochromatic  waves  1 29-34) . However,  it  is  ol 
importance  to  study  the  transient  response  ot  these  structures  to  finite  pulses.  In  tins  section 
the  reflection  and  transmission  of  a rectangular  pulse  and  a Gaussian  pulse  impinging  on  a 
periodic  structure  of  finite  length  C are  studied.  In  the  analytic  analysis,  the  general  case  ol 
transversely  bounded  structures  is  considered  wherein  the  coupled  waves  could  be  diflcrcnt 
modes  and  therefore  have  different  longitudinal  wave  vectors  [34) . For  simplicity,  and  with 
no  loss  of  generality,  the  numerical  examples  will  correspond  to  the  special  case  ol  trans- 
versely unbounded  structures.  The  eoupled-wave  approach  135-37)  is  used  and  its  limita- 
tions are  briefly  discussed  relative  to  the  more  general  Floquet  approach  (29-33).  There  is 
no  limitation  on  the  type  of  periodicity  as  long  as  the  coupl.ng  eoetticient  x is  small  enough 
to  satisfy  the  coupled-wave  approach  conditions. 


TRANSFER  FUNCTIONS 

Let  us  consider  a wave  F(z)  Up(x)e  P propagating  in  a transversely  bounded 
periodic  structure.  The  index  p represents  the  guided  mode,  0p  is  its  longitudinal  wave  vec- 
tor, and  Up(x)  is  the  transverse  dependence.  This  wave  can  be  eflectively  coupled  to  a 

backward  q mode  wave  B(z)  U(j(x)e  P it  the  phase-matching  condition, 

0p  + 0q  = 2»r/A.  (28 

is  satisfied,  where  A is  the  period  of  the  structure  (figure  47).  Let  be  the  frequency  at 
which  phase  matching  occurs.  The  coupled-wave  equations  are  well  known  to  be  1 35-37 1. 


and 
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where  Xpq  IS  the  coupling  coefficient  which  depends  on  the  transverse  I unctions  Up(x)  and 
U (x)  (34) , 6p  and  5t}  are  the  wave  vector  mismatch,  and  they  are  related  to  the  frequency 
match  by 
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Figure  47.  Cun  figuration. 
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where  Au>  = to  - , 

oj  = angular  frequency  , 

and  i//p  and  i//q  are  the  slopes  of  the  dispersion  diagram  of  the  guided  wave.  If  C is  the  length 
of  the  coupling  region,  then  the  longitudinal  boundary  conditions  are  (see  figure  47) 

R-  C/2)  = I , F(C/2)  = T , (31) 

B(-  C/2)  = R , and  B(G/2)  = 0 , (32) 

where  T and  R are  the  transmission  and  reflection  coefficients.  The  solution  for  the  coupled- 
wave  equation  (29)  with  the  above  boundary  condition  gives 

K r-T-  03) 

0 + 0fl 

DcothDC-  i 


and  T =- 


i,s«' V* 


5n  + 6u 

D cosh  DC  - i 4-^ — - sinhDC 


where 


If  the  two  coupled  modes  are  identical,  then  5p  - 6(j  - 5, 


(34) 


IX 

(35) 

DcothDC -ifi  ' 

D 

(36) 

1)  oshDC  - i5  sinDC 

i i 

and  D = (x“  - 

from  equations  (33  and  (34)  we  can  derive  the  chaiacteristie  of  the  periodic  structure: 


Maximum  reflection: 

Rm  = i tanh(XpqC) 

(37) 

Minimum  transmission: 

Tin  1 /cosh(xpqC) 

(38) 
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3. 


(39) 


R - 0 and  T - 1 for:  §p  + 6(j  = 2(x“  + (rur/C)-]  n = 1,2,...; 

the  corresponding  phase  of  R is  equal  to  (2n  - 1 ) 7r/2 

4.  |R|-  + |T|-  = 1 

In  figure  48  are  plotted  the  magnitudes  of  the  reflection  and  transmission  coefficients  as  a 
function  ot  5G  for  a given  value  of  xE  Figure  49  shows  the  equimagnitude  and  equiphase 
curves  ot  the  reflection  coefficient  plotted  as  a function  of  xC  and  6C.  Figure  48  corresponds 
to  a horizontal  line  (x^  = 2)  across  figure  49. 


TRANSIENT  RESPONSE 

The  reflection  and  transmission  of  a pulse  with  a relatively  narrow  spectrum  near 
uipq  can  now  he  derived  trom  the  reflection  and  transmission  coefficients  of  each  spectra 
component.  For  an  incident  pulse  g(t)  of  spectrum  G(co)  , 


G(co)  = 


♦ OO 

/g(t)e  iojt  dt 


(40) 


the  reflected  pulse  is 


+ 0 

f 


G(cj)  R(cj)c  iwt  dco 


(41) 


and  the  transmitted  pulse  is 


(42) 


Numerical  inversion  has  been  used  by  a number  of  authors  [38,  39  and  others]  to  obtain 
time  response  to  radiation,  propagation  and  scattering  problems  by  convolution  or  transforms. 
The  method  used  here  is  the  Cooley-Tukcy  fast  Fourier  transform  (FFT)  [40],  with  21 1 = 2048 
samples  in  the  FFT  to  calculate  equations  (41 ) and  (42). 

The  reflection  and  transmission  of  several  rectangular  pulses  were  studied,  with 
carriei  frequency  cjq  = ojpq  oj|,  and  where  toj  and  to  , correspond  to  the  first  two 
zeros  ol  the  reflection  coefficient.  The  normalized  pulse  lengths  r were  chosen  as  0.25  time 
unit  and  2.0  time  units,  where  each  time  unit  corresponds  to  the  transit  time  of  the  pulse 
across  the  slab  of  length  C. 

Figure  50  displays  the  reflected  and  transmitted  pulses  which  result  from  an  incident 
rectangular  pulse  of  r = 0.25  with  a carrier  frequency  exactly  at  phase  match  (to  = coq). 


REFLECTION  COEFFICIENT  AMPLITUDE  TRANSMISSION  COEFFICIENT  AMPLITUDE 


NORMALIZED  FREQUENCY  bV. 

F igure  48.  Magnitude  of  the  reflection  and  transmission  coefficients  as  a function  of  the 
normalized  frequency  mismatch:  5C=  AtoC/c  for  a fixed  coupling:  -■ 
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figure  49.  Equimagnttude  and  equiphasc  curves  of  the  reflection  coefficient  as  a function  of  6C 
and  X^-  The  reflection  curve  in  figure  48  corresponds  to  a horizontal  line  (x?  = 2)  across  the  grapl 
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Figure  50.  Reflected  and  transmitted  pulses  for  different  values  of  the  coupling  coefficient 
The  incident  rectangular  pulse  has  a length  of  0.25  time  unit  and  a carrier  frequency  cog  at 
the  Bragg  condition. 
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1 mcm"  Val"CS  °f  COl,plinB  X-  Nwe  ^ »a»  coupling 
nb)  tl,L  pu,SL  ,s  sprcad  ovcr  two  time  units.  This  is  to  be  expected  because  the 
uiergy  ts  reflected  from  successive  stratifications  and  the  echo  from  the  last  layer  would 

tr  ^ l*  r t V * * tF|iP  tllllC  °f  tW°  llll‘tS'  Sincc  t,le  reflections  arc  relatively  weak  (yC  « 1 ) the 
ransmttted  pulse  .s  similar  to  the  delayed  incident  pulse.  As  the  coupling  is ^incaSd  (stronger 
reflections  from  each  layer),  multiple  interference  plays  an  appreciable  role  which  leads  to 
wiH  |,Cn  S "I"  a lo,,gCT  time  P"iod  <41  >•  Simultaneously,  the  reflected  and  traTsm'^  pulses 
pull!  U«  r,’ k"’tl 'n  Al  hiel’  coupling  (figure  50,  j)  the  refleclcd 

n i„  . pu  j 'me  unit  and  decreases  rapidly  afterward  because  most  of  the 

/ ^ 'o'  bCCU  rC'lcCtcd  back  by  the  first  reflection.  However,  a relatively  complicated 

interference  tad  of  multiple  pulses  is  present.  The  transmitted  pulse  in  this  III  coZn  two 
main  narrow  peaks  which  correspond  to  differentiation  of  the  incident  pulse 

len,th  r - !hC  SamC  SL‘q"enCC  ,s  rcpcatcd  for  the  ^se  of  a rectangular  pulse  of 

re  si  7-S  ° L1;iracteristies  should  be  noted.  At  low  coupling  the  reflected  pulse 

(figure  51a,  ,,  similar  to  the  autocorrelation  of  the  incident  pulse.  This  is  due  to  the  sTm  ilarit  v 

the  tec,ion  cocffidcnt  in  thc  ^ rhry 

simil  ev  , , • ' 7 ‘n  8CnCral-  A‘  high  coliplille<  ‘he  transmitted  pulse  is  again 

m.lar  (except  tor  wave  interference  echoes)  to  the  derivative  of  the  incoming  pulse  As  in 

the  previous  case,  the  transmitted  pulse  is  similar  to  the  input  pulse  for  weak  coupling 

of  width?-"  I’2  dlSp,;!yS  thct  ;e,lectcd  aild  transmitted  pulses  for  a Gaussian  incident  pulse 

nentso!  th~  r " ^ ,sUakcn  Ut ,hc  l/e  vahies)‘  Since  the  high-frequency  compo- 

llc  ^'auss,an  pulse  ;,re  waller  than  those  of  the  rectangular  pulse  the  characteristic 
s larpen.ng  of  the  edges  of  the  transmitted  pulse  found  in  previous  cases  is  not  found  here 

HoweveMf  he' Cm,  !"  figUrc  52'  a11  0lltput  tend  to  mimic  the  input  pulse. 

ma  t f ' lan  pu,sc  ,s  narr°wed  to  r = 0.25  (not  shown  here),  the  results  are 

similar  to  a rectangular  pulse  of  the  same  width. 

I if,  . j i'ial,y’.  a S^ies  of  a,rvcs  was  generated  for  the  case  in  which  the  carrier  frequency  is 
narrow  (r  0 5?  54)  .cro  of  the  reflection  coefficient  Both 

‘ , "5  ‘“ld  Wldc  (T  ~ -0)  rectangular  pulses  were  considered.  In  these  cases  the 

rc  lectec1  output  is  greatly  distorted  by  a number  of  auxiliary  pulses  including  some  that  arc 

mg  Pll,SC-  AS  bCf°rC- thC  ,rallSmittL‘d  Pld-  generafshape  if  the 


LIMITATIONS  AND  APPLICATIONS 

the  main  fWmd^n  t'h  ' * ^ °"  ‘'1°  couP,cd-wave  aPProach,  which  only  accounts  for 

the  main  forward  and  backward  waves  and  neglects  the  presence  of  higher  harmonics  This  is 

ov.!r  a^)vnt! US  1C  COUpI,ng  IS  wcak  cnou8h  so  t're-  wave  amplitude  does  not  change  much 
a v a vc  length  or  over  a structure  period,  and  the  input  bandwidth  is  not  too  large  Other 

Z- 1! .T*  air,acl: 1 2M31  whidi  *»»•* Z lot  Z 

left  for  w7k  1 " ,SS",,|,|C  10  f0,ra,"a'c-  is  rcla,ivc|y  numerically  and  is 

Periodic  structures,  in  different  forms,  have  played  an  important  role  in  manv  fields- 

!!^iTn!k^7!!!i  ?;  and  SO"df  a'c  *«*«%.  *he«r  upplicaliot, in  the 

ergmg  field  of  opt, cal  communication  and  integrated  optics  [33]  has  attracted  appreciable 

■:tz  ^irv"ai  T”1  ton"m"’ica,io"  wmi  be  *»  dwa|  Cl? 

1 aper  arc  relevant.  For  instance,  a periodic  structure  can  be  used  as  a “derivator”  for  edge 
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1 


NORMALIZED  TIME 


f igure  52.  Reflected  and  transmitted  pulses  for  the  case  of  an  incident  Gaussian  pulse  of  width  2.0 
time  units  and  carrier  frequency  ojq. 
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REFLECTED  OUTPUT  TRANSMITTED  OUTPUT 


Figure  53.  Reflected  and  transmitted  pulses  for  the  case  of  an  incident  rectangular  pulse  of 
widths  0.25  time  unit  and  2.0  time  units  and  carrier  frequency  CO]  at  the  first  ze.'<  of  the 
reflection  coefficient. 
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detection  (figure  5 lj)  or  multiple-pulse  generator  (figure  50j).  The  above  results  can  also  be 
applied  to  such  fields  as  microwave  filters,  optical  filters,  and  sounders  of  subsurface  layers. 
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